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113  ABSTRACT  1 

A mathematical  model  of  an  engagement  between  a Navy  ship  and  a number  of 
anti-ship  cruise  missiles  is  described.  The  model  is  designed  to  be  used  in  Navy 
planning  activities  involving  ship/cruise  missile  engagements. 

The  engagement  model  described  portrays  the  ship's  defensive  weapons  as 
being  rapid-fire  guns  and/or  SAM  systems,  and  allows  for  the  portrayal  of  enemy 
jamming  or  other  ECM  activity.  User-specified  notions  of  partial  damage  to  the 
ship's  defensive  weapons  can  also  be  accommodated.  The  outputs  of  the  model  are 
probability  measures  that  the  ship  and  each  of  the  cruise  missiles  survive  as 
functions  of  time.  Inputs  to  the  model  are  detailed  performance  characteristics 
of  sensor  systems,  weapon  systems  and  cruise  missiles.  It  is  intended  that  these 
inputs  be  supplied  where  necessary  by  engineering  models,  one  of  which  is  described 
in  this  report.  The  model  is  analytic  in  nature  in  that  the  output  probability 
measures  are  obtained  by  numerically  solving  equations  portraying  the  relationships 
between  the  processes  which  make  up  the  engagement. 
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FOREWORD 


This  report  documents  research  conducted  for  the  Office  of  Naval 
Research  by  Vector  Research,  Incorporated*  ( VRI ) under  contract 
N00014-72-C-0300.  The  Research  is  a continuation  of  earlier  activities 
conducted  by  VRI  under |the  same  contract  and  described  in  the  report 

entitled  Development  of  Analytic  Methodoloqy  for  Naval  Planning  Areas 

i 

(VRI  report  number  0N$-1  FR  73-1).  One  result  of  the  early  research  was 

/ 

a finding  that  an  analytic  model  of  an  engagement  between  a ship  and  a 

( 

number  of  cruise  missiles  would  be  useful  as  an  aid  to  naval  planning. 
Such  a model  has  b?en  developed  and  is  described  in  this  report.  The 
model  is  mathematical  in  character.  This  report  is  written  for 
technical  personnel  who  will  program  it  and  assist  in  its  application. 
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' 1.0  INTRODUCTION 

The  research  documented  in  this  report  represents  a continuation  of 
work  begun  in  1972  which  had  as  objectives: 

(1)  the  study  of  Navy  planning  areas, 

(2)  the  examination  of  current  models,  and 

(3)  the  development  of  analytic  structures. 

Initial  work  in  pursuit  of  these  objectives  focused  on  airborne  anti- 
submarine warfare  and  on  fleet  air  defense;  models  utilized  by  the  Navy  in 
conducting  planning  in  both  areas  were  examined.  Two  results  of  this  exam- 
ination were  the  observations  that  the  most  commonly  used  models  were  Monte 
Carlo  simulations  and  that  there  was  a strong  possibility  that  analytic 
models  could  be  developed  as  supplements  to  or  substitutes  for  such  models. 
Following  discussions  with  personnel  at  ONR  and  the  Office  of  the  Chief  of 
Naval  Operations  (OP  96),  it  was  decided  to  concentrate  activities  in  the 
area  of  fleet  air  defense,  in  particular  fleet  air  defense  against  cruise 
missiles,  because  of  the  threat  posed  by  such  missiles  and  because  of  the 
fact  that  existing  models  were  of  only  limited  applicability  in  studying 
ship/cruise  missle  engagements. 

The  remainder  of  the  first  year's  efforts  were  concerned  with  the 
development  of  an  analytic  model  of  an  engagement  between  a single  ship 
and  a number  of  cruise  missiles.  In  addition,  preliminary  work  was  carried 
out  to  determine  the  structural  requirements  of  a fleet/multiple  cruise 

missile  engagement  model.  A hierarchical  analytic  structure  was  proposed. 

y 

i.  In  this  structure  one  set  of  models,  referred  to  as  engineering  models, 

are  used  to  map  subsystem  hardware  characteristics  into  subsystem  perform- 
ance characteristics.  A second  set  of  models,  referred  to  as  operational 
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models,  are  used  to  integrate  subsystem  performance  characteristics  to  pro- 
vide measures  of  total  system  effectiveness  (of  either  a ship  or  a fleet) 
in  an  engagement  with  multiple  cruise  missiles.  The  advantages  of  this 
approach  which  are  discussed  fully  in  [Bonder,  Cherry  and  Miller,  1973], 
include  not  only  economy  of  operation  but  also  flexibility  in  that  pro- 
posed hardware  changes  can  be  more  easily  evaluated  and  areas  in  which  im- 
provement is  desirable  or  necessary  can  be  more  easily  identified. 

Work  was  carried  out  during  the  first  year's  contract  on  both  opera- 
tional and  engineering  models.  A kill -rate  structure  was  adopted  for  the 
single  ship/multiple  cruise  missile  operational  model  and  rapid  fire  gun 
system  kill  rates  were  developed.  The  extension  of  such  a structure  to 
a fleet  operational  model  was  discussed  and  evaluated.  Requirements  for 
engineering  models  were  identified  and  modeling  was  carried  out  in  the  areas 
of  radar  power  management,  coherent  and  noncoherent  detection,  Doppler 
discrimination  and  target  tracking  and  prediction.  The  results  of  devel- 
opment activities  in  both  operational  and  engineering  models  are  described 
in  [Bonder,  Cherry  and  Miller,  1973]. 

Subsequent  to  the  completion  of  the  activities  of  the  first  year, 
initial  findings  and  results  were  reviewed,  and  the  decision  made  to  con- 
centrate the  second  year  of  the  research  activity  on  the  development  of 
analytic  operational  models  of  engagements  between  cruise  missiles  and 
single  ships.  A model  representing  engagements  between  a single  ship 
and  a number  of  cruise  missiles  was  developed  as  a result  of  that  effort 
and  is  described  in  this  report. 

The  single  ship/multiple  cruise  missile  engagement  model  is  an  1 

c 

operational  model  (i.e.,  it  maps  performance  descriptions  for  sub- 
systems into  measures  of  effectiveness  of  the  engagement)  and  ■ 
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is  analytic  in  character.  The  principal  mathematical  structures  are  sto- 
chastic, and  the  outputs  are  probabilistic  descriptions  of  measures  of 
effectiveness  (outcomes  of  engagements)  as  functions  of  time.  The  elements 
and  processes  represented  in  the  model  include: 

• the  ship's  defensive  weapon  systems  and  weapon  assignment  logic, 

• damage  to  defensive  weapons, 

• in-flight  destruction  of  cruise  missiles, 

• jammers  and  other  ECM  systems  employed  against  defensive  weapon 
systems,  and 

• raid  characteristics. 

The  engagement  model  requires  as  inputs  the  performance  characteristics 
of  systems  (missiles,  defensive  weapons,  etc.)  including  such  items  as  le- 
thality characteristics  of  cruise  missiles  with  respect  to  the  ship  (and 
vice  versa),  tracking  error  characteristics  of  the  ship's  radar  in  an  EW 
environment,  etc.  While  it  is  anticipated  that  some  of  the  model's  inputs 
will  be  determined  directly  from  existing  data,  it  is  believed  that  other 
inputs  should  be  obtained  as  outputs  of  engineering  submodels  that  map 
physical  characteristics  of  subsystems  into  performance  descriptors  of 
those  subsystems,  which  can  then  be  input  to  the  engagement  model.1  Accord- 
ingly, in  addition  to  the  engagement  model  itself,  an  engineering  submodel 
that  permits  tracking  error  characteristics  of  radars  and  other  sensors  in 
an  EW  environment  to  be  determined  as  a function  of  the  physical  character- 
istics of  the  radar  was  also  developed  during  the  second  year's  effort  and 
is  documented  in  this  report. 

1 Wi th  this  arrangement  the  effect  of  variations  of  these  subsystems  on  the 
outcome  of  the  engagement  can  be  examined  by  varying  the  inputs  to  the 
submodels  --  a convenient  and  efficient  approach  which  may  facilitate  the 
search  for  improved  systems. 


The  remainder  of  the  report  consists  of  two  chapters  and  three 
appendices.  The  single  ship/multiple  cruise  missile  engagement  model 
is  described  in  chapter  2.0.  Expressions  for  the  kill  rates  and  weapon 
firing  probabilities  (both  of  which  are  required  for  the  overall  model) 
in  terms  of  subsystems  performance  descriptors  aie  also  developed  in 
chapter  2.0.  Suggestions  for  refining  the  model  and  recommendations 
for  extending  the  model  to  represent  multiple  ship/multiple  cruise 
missile  engagements  are  discussed  in  chapter  3.0.  An  engineering 
model,  which  may  be  used  to  generate  some  of  the  parameters  involved 
in  the  expressions’ for  kill  rates,  is  described  in  appendices  A and  B. 

A system  of  differential  equations  employed  in  the  overall  model  is 
analyzed  in  appendix  C. 
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2.0  AN  ANALYTIC  MODEL  OF  A SINGLE  SHIP/MULTIPLE 
CRUISE  MISSILE  ENGAGEMENT 


An  analytic  model  of  a single  ship/multiple  cruise  missile  engage- 
ment is  described  in  detail  in  section  2.1,  together  with  the  principal 
assumptions  used' in  the  development  of  the  model.  Expressions  for  "kill 
rates"  and  "firing  probabilities,"  which  are  fundamental  to  the  model, 
are  developed  in  sections  2.2  and  2.3  respectively. 

2. 1 Overall  Structure  of  the  Model 

2.1.1  Principal  Assumptions  and  Inputs, -to.  the  Model 

A single  ship/multiple  cruise  missile  engagement  is  viewed  as 
taking  place  in  two.  phas.es.;  an  ea*Ty  phase  in  which  the  cruise  missiles 
are  engaged  only  by-  defans iv.O'^aif craft  at  a disla.ice  from  the  ship,  and 
a late  phase  in  whichotite  cruise  missiles  are  engaged  only  by  the  ship's 
defensive  weapons at  relatively  close  ranges.  Because  defensive  air- 

S' 

craft  wiU- seldom  be  used  at  the  same  time  as  the  ship's  defensive 

s' 

weapons,  the  view  is.  that  no  such  aircraft  participate  in  the  late 
phase  of  the  engagement. 

The  engagement  model  described  in  this  report  portrays  only  the 
late  phase  of  the  engagement1  of  which  figure  1 is  an  illustration. 

^ther  models  which  could  portray  the  early  phase  of  the  engagement 
(aircraft  vs  cruise  missiles)  are  available  elsewhere. 


* O' 
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The  following  principal  assumptions  about  this  phase  have  been  made:1 

(1)  The  ship  detects  cruise  missiles  at  deterministic  ranges  from 
the  ship.2  Detection  ranges  may  vary  with  the  direction  of 
incidence  of  the  cruise  missiles.  Upon  detecting  a cruise 
missile,  the  ship  initiates  tracking  and  continually  updates 
its  estimates  of  what  the  position  of  the  cruise  missile  will 
be  at  future  times. 

(2)  The  ship's  defensive  weapons  are  limited  to  guns  and  SAM  sys- 
tems. 

(3)  The  ship  can  redirect  in-flight  rounds  from  both  guns  and 
SAM  systems  at  any  time  prior  to  their  arrival  at  the  cruise 
missile  at  which  they  have  last  been  targeted. 

(4)  The  ship  always  assigns  all  available  defensive  weapons  to 
that  surviving  cruise  missile  which  will,  if  not  sooner  de- 
stroyed, be  the  next  to  impact  the  ship. 

The  above-listed  principal  assumptions  serve  to  bound  (and,  in  the  case 
of  assumption  (3),  simplify)  the  scenario  to  be  modeled  as  well  as  to 
define  some  of  the  types  of  input  data  which  will  be  required  by  the 
model.  The  manner  in  which  each  assumption  accomplishes  these  objec- 
tives is  briefly  discussed  below. 


^ther  assumptions,  more  minor  in  character  than  those  listed  here,  will 
be  introduced  and  used  as  needed. 

2This  assumption,  and  the  others  as  well,  can  be  relaxed.  Methods  for 
relaxing  the  assumptions  will  be  discussed  later. 


: i.  .r.-WMwa&. 
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Detection  and  Tracking.  Since  the  ranges  at  which  a ship  can  detect 
incoming  cruise  missiles  can  be  predicted  with  reasonable  confidence,1 
these  ranges  have  been  viewed  as  being  deterministic  in  character.  The 
range  at  which  detection  of  each  cruise  missile  occurs,  together  with  the 
actual  flight  paths  and  velocities  of  the  cruise  missiles,  are  required 
as  input  to  the  model  in  order  to  obtain  the  time  which  the  ship  has  to 
engage  and  destroy  each  missile  before  it  impacts  the  ship  and  to  fix  the 
geometry  of  the  engagement. 

In  practice,  it  will  be  convenient  to  represent  the  actual  flight 
path  of  each  cruise  missile  by  a sequence  of  straight  line  segments  and 
constant-radius  turns  of  a specified  number  of  degrees.2  The  speed  of 
the  cruise  missile  is  taken  to  be  constant  on  each  segment. 

Ship's  Defensive  Weapon  Systems.  A single  ship's  present  defenses 
against  cruise  missiles  consist  of  rapid-fire  guns  and  SAM  systems.  The 
number  of  defensive  weapon  systems,  along  with  certain  performance  data 
which  affect  their  lethality  characteristics  with  respect  to  cruise  mis- 
siles, is  required  as  input  to  the  model. 


Perhaps  the  simplest  method  to  obtain  the  detection  range  is  to  use  the 
results  of  tests  of  the  ship's  missile  detection  radar  against  cruise- 
missile-like  objects  under  conditions  of  the  type  a user  desires  to  rep- 
resent in  the  model  (e.g.,  presence  or  absence  of  enemy  jamming,  weather 
conditions,  the  altitude  of  the  incoming  cruise  missiles,  etc.).  Alter- 
natively, the  ship's  detection  range  could  be  obtained  from  analytic 
models  (see,  for  example,  [Brennan  and  Hill,  1964]  and  [Kirkwood,  1965]). 
Similarly,  information  about  flight-path  estimation  for  incoming  cruise 
missiles  could  be  obtained  by  testing  the  equipment  involved  or  by  ana- 
lytic modeling.  An  approach  of  the  latter  type  is  documented  in  [Bonder, 
Cherry  and  Miller,  1973]. 


Procedures  of  this  type  for  handling  flight  paths  have  been  developed  at 
the  Systems  Research  Laboratory  at  the  University  of  Michigan;  they  are 
documented  in  [SRL,  1969]. 


■ i n iiVi*  i 


Fire  Redirection.  On  the  one  hand,  it  is  sometimes  possible  for  a 


ship  to  redirect  an  in-flight  SAM  from  the  cruise  missile  at  which  it  was 
initially  targeted  to  another  cruise  missile.1  On  the  other  hand,  redirec- 
tion of  in-flight  gunfire  by  the  ship  does  not  appear  to  be  feasible. 
Nevertheless  it  is  here  assumed  for  mathematical  convenience  and  simplicity 
that  redirection  of  both  gunfire  and  SAMs  is  possible.  Strictly  speaking, 
therefore,  this  assumption  has  the  effect  of  making  the  engagement  model 
described  here  generate  an  upper  bound  on  the  effectiveness  with  which 
real  ships  can  fight  cruise  missile  engagements.  While  an  analysis  of 
such  engagements  is  possible  without  this  assumption,  it  is  not  clear  that 
the  additional  realism  thus  gained  would  be  worth  the  concomitant  increase 
in  mathematical  and  computational  complexity. 

In  view  of  the  assumption  that  the  ship  has  the  capability  to  redirect 
in-flight  fire  from  any  weapon  system,  it  is  convenient  to  say  that  a weap- 
on system  "is  firing  at"  a cruise  missile  at  a particular  time  whenever 
fire  from  that  weapon  is  arriving  at  the  cruise  missile  at  that  time. 

Ship's  Weapon  Assignment  Logic.2  Although  other  weapon  assignment 
logics  are  available,  the  logic  described  in  assumption  (4)  above  closely 
approximates  one  which  has  been  used  in  practice.  This  logic  has  there- 
fore been  assumed  as  an  initial  basis  for  modeling. 

Redirection  of  a SAM  is  usually  possible  in  practice  until  the  SAM  enters 
a terminal  lock-on  phase. 

2Ship's  weapon  assignment  logics  in  general,  and  the  one  assumed  for  the 
engagement  model  in  particular,  are  discussed  at  greater  length  in 
section  2.3. 
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Notice  finally  that  some  of  the  above  assumptions  state  that  certain 
aspects  of  the  engagement  are  regarded  as  being  deterministic  instead  of 
stochastic.  For  example,  assumption  (1)  says  that  the  ranges  at  which  the 
ship  detects  cruise  missiles  is  here  taken  to  be  deterministic,  whereas 
these  ranges  are  often  treated  stochastically  elsewhere  (i.e. , in  ether 
models).  These  assumptions  may  therefore  be  seen  as  conditioning  the  en- 
gagement upon  the  indicated  types  of  information  being  fixed.  This  con- 
ditioning could  be  removed,  if  desired,  by  embedding  the  engagement  model 
to  be  described  here  in  a richer  structure  in  which  these  deterministic 
aspects  would  instead  be  treated  stochastically.  This  would  entail  an 
increase  in  mathematical  and  computational  complexity.  Since  much  of  the 
complexity  of  the  interactions  between  the  ship  and  the  cruise  missile  can 
already  be  analyzed  under  the  assumptions  listed  above,  it  is  desirable 
to  postpone  enrichments  of  the  types  suggested  above  until  such  time  as 
experience  gained  with  the  existing  model  indicates  that  they  are  needed. 

2.1.2  Outputs  and  Overall  Mathematical  Structure  of  the  Model 

The  outputs  of  the  single  ship/multiple  cruise  missile  engagement 
model  are  of  two  types: 

• the  conditional  probability  that  each  cruise  missile  survives  the 
engagement  up  to  time  t,  given  that  the  ship  has  sustained  a cer- 
tain degree  of  damage  at  that  time  (the  output  is  provided  for 
each  time  t and  each  possible  degree  of  damage  to  the  ship),  and 

• the  probability  that  the  ship  sustains  each  possible  degree  of 


damage  at  each  time  t. 
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The  overall  structure  of  the  engagement  model  is  based  on  renewal 
theory.  Fundamental  to  the  model  are  "kill  rates"  and  what  are  here 
called  "firing  probabilities."  The  kill  rate  associated  with  each  defen- 
sive weapon  system/cruise  missile  pair  at  time  t is  the  conditional  prob- 
ability that  the  weapon  system  will  kill  the  cruise  missile  in  a small 
time  interval1  [t,  t+6),  given  that  the  cruise  missile  is  still  alive  at 
time  t and  that  tjie  weapon  is  firing  at  the  cruise  missile  at  that  time.2 
The  firing  probability  p(t)  associated  with  each  defensive  weapon  system/ 
cruise  missile  pair  is  the  probability  that  the  weapon  is  actually  firing 
on  the  missile  at  time  t.3 

It  is  convenient  to  number  the  cruise  missiles  in  the  order  in 
which  they  will,  if  not  sooner  destroyed,  impact  the  ship.  Thus  cruise 
missile  1 is  that  missile  which,  if  not  destroyed,  will  impact  the  ship 
before  any  of  the  other  cruise  missiles;  cruise  missile  2 is  that  missile 
which  will,  if  not  destroyed,  be  the  next  soonest  to  impact  the  ship, 
and  so  on.  The  convention  here  is  that  time  is  measured  from  zero 
starting  with  the  time  at  which  the  cruise  missile  1 is  detected  by  the 
ship.  The  time  at  which  cruise  missile  j will,  if  not  sooner  destroyed, 
impact  the  ship  is  denoted  by  "t^,"  so  that  0 _<  t-j  £ t£  £ follows 

from  these  conventions.  The  ship's  defensive  weapon  systems  may  be  numbered 
in  an  arbitrary  order. 


iHere,  as  elsewhere  in  this  report,  5>0  is  a small  time  interval. 

2Recall  that  a weapon  is  said  to  be  firing  at  a cruise  missile  at  time  t 
if  and  only  if  lethality  from  the  weapon  is  arriving  at  the  missile  at 
that  time. 

3The  reason  why  the  firing  of  a defensive  weapon  upon  a particular  cruise 
missile  is  treated  stochastically  will  be  explained  later. 


I 


Unlike  the  cruise  missiles,  the  ship's  defensive  weapons  are  viewed 
as  being  capable  of  sustaining  partial  damage  (as  would  happen,  for  example, 
when  guns  switch  to  optical  tracking  because  of  damage  to  the  radar  system). 
Partial  damage  to  the  ship's  defensive  weapons  can  conveniently  be  de- 
scribed by  introducing  and  using  "damage  categories"  for  the  weapons  and 
"damage  states"  for  the  ship.  The  damage  categories  for  one  of  the  ship's 
defensive  weapons  are  statements  that  are  descriptive  of  that  weapon's 
capability  to  continue  to  function  in  the  engagement  and  may  be  chosen 
arbitrarily  by  the  user  of  the  model.  For  example,  the  damage  categories 
for  a weapon  might  be  chosen  as  follows1: 


damage 

category  number 


weapon  capability 
none  (weapon  destroyed) 
marginal 
medium 
almost  full 
full  (no  damage) 


The  damage  categories  for  the  other  defensive  weapons  may  be  chosen  simi- 
larly. If  there  are  n^  damage  categories  associated  with  defensive  weapon 
i,  and  if  there  are  W defensive  weapons  on  the  ship,  then  there  are  JJn.j 
combinations  of  damage  categories  that  could  characterize  the  status^f 
the  ship's  defensive  weapons  at  any  time.  It  is  convenient  to  call  each 
such  combination  a damage  state  of  the  ship.  Damage  states  are  denoted 


^he  terms  which  define  the  categories  must,  of  course,  themselves  be 
defined. 
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by  vectors  d,  where  the  i^  component  of  a damage  state  vector  d is 
the  number  of  the  damage  category  of  defensive  weapon  .ystem  i.  The  par- 
ticular damage  state  corresponding  to  all  of  the  ship's  defensive  weapons 
being  functional  at  full  capability  is  denoted  by  u.  At  time  0,  the  ship 
is  assumed  to  be  in  damage  state  u. 

The  notation  used  is  summarized  below: 

M = number  of  cruise  missiles  in  the  engagement; 

W = number  of  defensive  weapon  systems  on  the  ship; 
t.  = time  at  which  cruise  missile  j will,  if  not  sooner 

J 

destroyed,  impact  the  ship  (1  £ j £M); 
m j ( t | d ) = conditional  probability  that  cruise  missile  j will 

be  alive  at  time  t,  given  that  the  ship  is  in  damage 
state  d at  that  time  (1  £ j £ M,  all  d); 

•S'k-jjUld)  = conditional  probability  that  defensive  weapon  system 
i will  kill  cruise  missile  j in  the  interval  [t,  t+s), 
given  that  cruise  missile  j is  still  alive  at  time  t, 
that  the  ship  is  in  damage  state  d at  time  t and  that 
defensive  weapon  i is  firing  at  cruise  missile  j at 
time  t (a  function  of  t for  each  1 £ i £ W,  1 £ j £ M 
and  each  d); 

Pij(t|d)  = conditional  probability  that  defensive  weapon  system 
i will  be  firing  at  cruise  missile  j at  time  t,  given 
that  cruise  missile  j is  still  alive  at  time  t and 
that  the  ship  is  in  damage  state  d at  that  time  (a 
function  of  t for  each  1 £ i £ W»  1 £ j £ M and  each  d); 
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i 


( t)  = probability  that  ship  will  be  in  damage  state  d at 
time  t (a  function  of  t for  each  d);1 
&ed  = conditional  probability  that,  given  that  the  ship 
is  in  damage  state  e,  it  will  be  in  damage  state  d 
after  the  impact  of  one  more  cruise  missile2  (all  e. 


d). 


In  addition,  for  each  j = 1 , 


, M,  let  t"  and  t.  denote  times  satis- 
J J 


fying  tT  < tj  < tt  which  are  both  very  close  to  t^.  Put  t~+1  = tM+^  = 
■ - • 


The  kill  rates  k..(t|d)  defined  above  express  the  lethality  charac- 
^ J 

teristics  of  defensive  weapon  system  i with  respect  to  cruise  missile  j 
at  time  t when  weapon  system  i is  firing  at  that  cruise  missile  at  that 
time.  Whether  or  not  a weapon  system  is  firing  at  cruise  missile  j at 
some  time  depends  (according  to  the  assumed  ship's  weapon  assignment 
logic)  on  whether  or  not  cruise  missiles  1,  2,  •••,  j-1  have  been  de- 
stroyed or  have  impacted  the  ship  by  that  time.  Therefore,  since  cruise 
missile  survival  is  here  treated  probabilistically,  so  too  must  be  the 
assignment  of  defensive  weapons  to  cruise  missiles  at  any  time.  The  fir- 
ing probabilities  p..(t)  defined  above,  express  the  probability  that 

* J 


^he  damage  state  probabilities  q^-)  may  be  interpreted  as  ship  survival 
probabilities  as  soon  as  a suitable  subset  of  the  damage  states  is  identi- 
fied with  ship  survival  (and  the  complementary  subset  with  ship  non- 
survival). We  do  not  here  make  this  association,  preferring  instead  to 
leave  this  to  the  user  of  the  model.  The  qc|( • ) are  therefore  made  avail- 
able as  outputs  of  the  model  and  this  is  regarded  as  being  equivalent  to 
outputting  a ship  survival  probability  as  a function  of  time. 

2It  is  assumed  that  the  a ^ are  known  and  available  as  input  data. 
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in  terms  of  the  m.(t|d)  are  available.1  Specifically,  one  could  use 
J 

(3)  to  first  determine  the  m.(6|u)  for  1 < j < M.  This  information  can 

J 

then  be  used  to  determine  p. -(6 Ju)  for  1 < i < W and  1 < j < M.  This 

* J 

latter  information  could  then  be  used  in  (3)  to  obtain  the  ( 26 | u ) , 
etc.  Continuing  with  this  "bootstrap  technique"  would  yield  values  of 
the  m.(t|u)  and  p..(t|u)  at  any  desired  number  of  points  in  the  range 
0 < t < t-j . The  following  diagram  illustrates  the  order  of  the  computa- 
tions: 


m.(0|u)  m.(6|u)  m.(26|u) 


/ 


m.(tj|u) 


pij(°lu)  P-Ulu)  Pij(2«|u) 


pij(tT i u) 


The  analysis  just  given,  as  has  been  indicated,  suffices  to  deter- 
mine the  m.(t|d)  and  q,(t)  only  for  times  t < t-j . We  now  show  how  to 
compute  these  functions  for  times  t > t-j . To  this  end,  choose 
k(l  < k 5 M)  and  suppose  that  the  following  functions  and  quantities 
have  all  been  determined  as  indicated:2 


1 If  the  manner  of  dependence  of  the  firing  probabilities  p^ j ( 1 1 d)  on  the 
cruise  missile  survival  probabilities  m j ( t | d ) were  simple,  it  would  be 
possible  in  principle  to  insert  these  functions  in  (3)  (which  is  equiva- 
lent to  a system  of  differential  equations)  and  integrate  to  obtain 
formulas  for  the  m j ( • | d ) . An  approach  of  this  type  is  discussed  in 
appendix  C.  One  may,  however,  anticipate  that  the  integration  may  be 
difficult  if  the  p^j(t|d)  are  complicated  functions  of  the  m j ( • ! d ) and 
of  the  assumed  weapon  assignment  logic  for  the  ship.  In  view  of  this, 
it  is  better  to  plan  to  work  with  (3)  in  the  form  given  by  using  suit- 
able numerical  techniques. 

2This  assumption  is  true  when  k = 1 as  shown  above  (see  equations  (1) 
and  (2)  and  recall  the  bootstrap  technique  applied  to  (3)). 
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(a) 

qd(t)  for  all  d and  all  t <tfc; 

(b) 

qd(tp 

for  all  d; 

(c) 

nij{t|d) 

for  1 < j < M,  and  d and 

(d) 

mj-ft'ld) 

for  1 < j < m and  all  d; 

(e) 

Pij(t|d) 

for  1 < i < W>  1 < j < M 

The  probability  qd(tk)  that  the  ship  is  in  damage  state  d at  time  t+ 
is  given  jy  k 

qd(tk>  ' l qe(9  V^le)  ted  (4) 

for  all  d.  In  fact,  this  relation  persists  for  all  t satisfying  t*  < t < 


ck+l ' 


= z qe(t-)  mk(t-|e)  ted  ' (5) 

for  all  d and  all  t satisfying  t*  < t < tk+,.  In  particular. 


’d^M1  = l qe(tk)  «ed 


(6) 


for  all  d.  Notice  that  equations  (4),  (5),  and  (6)  are  all  computable  in 
the  sense  that  the  left-hand  side  of  each  can  be  determined  from  the 
corresponding  right-hand  side  because  all  of  the  right-hand  sides  involve 
only  probabilities  which  are,  by  the  above  assumption,  known.  Relation 
(5)  therefore  shows  that  part  (a)  of  this  assumption  remains  true  when  k 

is  replaced  by  M.  Similarly,  (6)  shows  that  part  (b)  of  the  assumption 
remains  true  when  k is  replaced  by  k+1. 

Next,  note  that  the  conditional  joint  probability  that: 

•the  ship  is  impacted  at  time  tfc,  and 


f 
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•the  ship  thereupon  goes  into  damage  state  d,  and 
•cruise  missile  j(k+l  i j i M)  survives  until  time  t£  ss  tj^, 
given  that  the  ship  is  in  damage  state  e at  time  tj^,  is  m^t^le)  • Aed. 


nj<tkld>  • <’d(tk>  = | "k^kl6)  2ed  qe(tk) 

+qd(t£)  • [mj(tjj|d)  - mk(tjj|d)] 


for  j = k+l,#",  M and  all  d,  whereas 

mj(tkld)  = 0 (») 

for  j=l,‘*‘,  k and  all  d because  cruise  missiles  !,•••,  k will  have  either 
been  destroyed  or  impacted  the  ship  by  time  t£.  Since  the  qd(t£)  which 
occurs  on  the  left  in  (7)  are  given  in  (4)  in  terms  of  probabilities 
evaluated  at  time  t7,  and  since  these  latter  probabilities  (together  with 
those  that  appear  on  the  right  in  (7))  are  all  known,  the  nij(t£|d)  on  the 
left  in  (7)  and  (8)  are  all  computable. 

Equation  (8)  persists  for  times  t satisfying  t£  £ t < t^+-j  because 
no  cruise  missiles  can  return  to  life  once  they  have  been  destroyed  or 
impacted  the  ship: 

m j ( 1 1 d ) = 0 (9] 

for  j=l , •••,  k,  all  d and  all  t satisfying  t£  <_  t < t^.  In  particular. 


I d)  • 0 


(10) 


c 
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for  j=l,  k and  all  d.  Furthermore,  an  argument  similar  in  every 
respect  to  that  given  before  shows  that 


i. 


ki- 


w 

m-(t+6|d)  = m-(tjd)  • [1  - 6 • l P^i(t|d)  k,,(t|d)]  (11) 

holds  for  all  j = k+1,  •••,  M,  all  d and  all  t satisfying  t^  <.  t < tfc+1. 

This  equation  can  be  numerically  bootstrapped  to  compute  the  cruise  missile 
survival  probabilities  m.(t|d)  and  the  weapon  system  firing  probabilities 

J 

Pi  j ( 1 1 d)  f°r  3^  t in  the  range  t£  <_  t < t^.  This  bootstrapping  begins 
with  the  choice  t = t^,  the  values  of  the  ( t^ | d)  which  will  then  be 

required  in  (11)  being  given  in  (7).  The  following  diagram  illustrates 
the  order  of  the  computations: 

mj<Wd) 

I 

»ij(tk+lld) 

In  particular,  the  m^- ( tj^+-j  | d)  for  i = k+1,  •••,  M and  all  d are  produced 
as  a result  of  this  procedure.  Thus,  equations  (9),  (10)  and  (11)  show 
that  parts  (c),  (d)  and  (e)  of  the  above  assumption  remain  true  when  k 
is  replaced  by  k+1. 

The  foregoing  discussion,  together  with  the  fact  that  the  above 
assumption  is  true  for  k = 1 , shows  by  induction  that  this  assumption 
is  true  for  all  k = 1,  •••,  M.  This  is  to  say  that  the  output  functions 


q,  (•)  and  m - ( * j d)  are  computable  for  all  damage  states  d and  all  cruise 
missiles  1 < j < H.  The  formulas  for  performing  these  computations  are 
those  which  have  been  developed  above.  The  procedure  may  be  recapped  as 
follows1: 

Step  (1).  Put  k = 1 and  compute 


qd(‘)  - 


1 if  d = u 
0 if  otherwise 


all  0 <_  t < t-j 


(12) 


qd(ti> 


1 if  d = u 
0 otherwise 


(13) 


Compute  also2  the  m.(t|u)  for  1 < j < H and  0 £ t < t,  by  numerically 

J * 

bootstrapping  the  following  differential  equation: 


m.(t+6|u)  = m^tju)  • [1  - 6 -4^pij.(t|u)  k..(t|u)]  . (14) 


i=l 


U 


Obtain,  in  particular,  the  values  m.(tT|u)  for  1 < j < M. 

J ' 

Step  (2).  Compute 

qd(tj)  = Lqe(tk)mk(tk|d)  £ed  a11  d’ 


(15) 


qd(t)  =Eqe(tjJ)mk(tjJ|e)  &ed  all  d,  all  tk  < t <tk+1. 


(16) 


Notice  that  the  computations  in  step  (2)  of  the  algorithm  should  be  per- 
formed in  the  order  indicated  because  the  qd(t£)  which  occurs  in  (18)  must 
first  be  computed  as  in  (15)  and  the  mi(t£|a)  which  will  be  required  to 
start  the  bootstrapping  in  (21)  must  first  be  obtained  as  in  (18). 


2Note  that  the  m ^ ( 1 1 d ) for  1 £ j £ M,  all  0 £ t < t-j  and  d 1 u are  undefined. 


£ 
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qd(tk+l)  = £<le(tk)mk(tk  e)  eed  3,1  d- 


(17) 


— V-  t&fcCt^le)  »ed  qe(t-) 


m.(t^|d)  -/  +<ld(tk)(mj(tkld>  - ",k<tkld)j;1 

all  k+1  < j < M,  all  d. 


(18) 


0 all  1 < j < k,  all  d, 
m j ( 1 1 d ) = 0 all  1 < j <_  k,  all  d,  all  t£  < t £ t^ , (19) 

mj(tk+l|d)  = 0 a11  1 iJik>  a11  d-  (20) 


In  addition,  compute  the  mj(t|d)  for  k+1  <_  j <_  M,  all  d and  all  t in  the 


range  t^  <_  t £ t^+-j  by  numerically  bootstrapping  the  following  differ- 
ential equation: 


m 


,(t+6|d)  =m,(t|d)  •[!-«•£  Pjjttld)  kj,.(t|d)].  (21) 


i=l 


ij 


Obtaii 


in,  in  particular,  the  values  (t^+1 |d)  for  k+1  < j < M and  all  d. 
Step  (3).  Increment  k by  1.  If  k >_  M,  stop;  otherwise  go  back  to  step 
(2). 

The  above  algorithm  constitutes  the  overall  mathematical  structure 
of  the  single  ship/multiple  cruise  missile  engagement  model.  In  order 
to  carry  out  the  above  procedure,  however,  it  is  necessary  to  first  deter- 


mine values  for  the  kill  rates  k - . ( 1 1 d ) and  the  firing  probabilities 

* J 


Pi j ( t I d) . Although  this  could  probably  be  done  experimentally,  it  would 


be  difficult  and  costly  to  do  so  and  would,  moreover,  provide  only  limited 
insight  as  to  possible  improvements  in  weapon  design  and  employment  tactics 


for  the  ship.  These  problems  can  be  avoided  by  using  analytic  models  to 
predict  the  kill  rates  and  firing  probabilities  in  terms  of  parameters 
which  characterize  the  cruise  missiles  and  the  ship's  defensive  weapon 
systems.  This  will  facilitate  using  the  overall  model  to  examine  promis- 
ing design  and  tactical  improvements.  Models  for  predicting  the  kill 
rates  and  firing  probabilities  have  therefore  been  developed  and  are  de- 
scribed in  the  next  two  sections. 

2.2  Determination  of  Kill  Rates 

The  overall  mathematical  structure  of  the  single  ship/multiple  cruise 
missile  engagement  model  described  in  the  previous  section  involved  certain 
"kill  rates"  which  characterize  the  capabilities  of  the  ship's  defensive 
weapon  systems  with  respect  to  cruise  missiles.  Each  kill  rate  k..(t|d) 

* J 

was  the  instantaneous  conditional  probability  that  defensive  weapon  system 
i will  kill  cruise  missile  j in  a short  time  interval  [t,t+6),  given  that 
cruise  missile  j is  still  alive  at  time  t,  that  the  ship  is  in  damage  state 
d at  time  t,  and  that  weapon  system  i is  firing  at  cruise  missile  j at 
time  t.  It  was  mentioned  that  kill  rates  for  both  types  of  the  ship's  de- 
fensive weapon  systems  — gun  emplacements  and  SAM  systems  — may  be  viewed 
as  being  functions  of  certain  parameters  associated  with  rhe  weapon  systems 
and  the  cruise  missiles.  The  dependence  of  kill  rates  for  rapid-fire  gun 
systems  on  such  parameters  has  been  described  by  [Bonder,  Cherry  and  Miller 
1973];  for  completeness,  these  ideas  are  discussed  in  section  2.2.1  below. 
Expressions  for  the  kill  rates  for  SAM  systems  as  functions  of  such  param- 
eters are  developed  in  section  2.2.2. 
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2.2.1  Kill  Rates  for  Rapid  Fire  Gun  Systems 

The  approach  used  here  to  determine  the  functional  form  of  the  kill 
rates  for  rapid-fire  gun  systems  involves  decomposing  the  process  by  which 
such  guns  kill  cruise  missiles  into  several  parts. 

Implicit  in  the  notion  of  a kill  rate1  k(t)  for  a particular  gun- 
system/cruise-missile  pair  at  time  t is  the  assumption  that  the  killing 
of  the  missile  by  rounds  from  the  gun  is  probabilistic  in  character.  If 
we  further  assume  that  the  killing  of  the  cruise  missile  at  any  time  t 
by  each  of  the  rounds  from  the  gun  which  arrive  at  the  cruise  missile  at 
time  t are  independent  events,  then  k(t)  may  be  written 

k(t)  = r(t) -SSKP(t) 


where: 

r(t)  = rate  at  which  rounds  from  the  gun  arrive  at  the 
cruise  missile  at  time  t;  and 
SSKP(t)  = single  shot  kill  probability  at  time  t. 

Note  that  the  rate  r(t)  at  which  rounds  arrive  at  the  cruise  missile 
at  time  t is  not  the  same  as  the  rate  at  which  the  rounds  were  fired  at  an 
earlier  time2  t-t^  (so  as  to  reach  the  missile  at  time  t)  because  of  Doppler 
effects  arising  from  the  motion  of  the  missile  relative  to  the  gun.  How- 
ever, since  the  rounds  from  the  gun  travel  so  much  faster  than  the  cruise 
missile,  these  Doppler  effects  will  be  small  and  may  therefore  be  neglected. 
The  above  equation  for  k(t)  then  becomes 

k(t)  = F-SSKP(t)  (22) 

[For  purposes  of  this  and  the  next  section,  it  is  convenient  to  replace  the 
notation  kjj(t|d)  for  a kill  rate  by  simply  k(t).  In  doing  this,  it  is 
understood  that  we  have  in  mind  a particular  defensive  weapon  system  i,  a 
particular  cruise  missile  j,  and  a particular  damage  state  d for  the  ship. 

2Here,  tf>0  represents  the  time  of  flight  for  rapid-fire  gun  rounds  arriv- 
ing at  the  cruise  missile  at  time  t. 


F = gun  system  firing  rate  ( a constant  which  appropriately  aver 

ages  the  burst  and  inter-burst  periods  for  the  gun);  and 

SSKP(t)  = single-shot  kill  probability  at  time  t. 

The  firing  rate  F is  here  taken  to  be  a basic  descriptor  of  the  gun 
system.  The  remainder  of  this  section  is,  therefore,  devoted  to  the  deter- 
mination of  the  single-shot  kill  probability  SSKP(t)  at  time  t. 

The  method  suggested  here  for  computing  single-shot  kill  probabilities 
for  rapid-fire  gun  systems  is  based  on  the  following  assumptions: 

• for  purposes  of  the  computation  of  the  SSKP(*)»  the  cruise  missile 
is  equivalent  to  an  effective  vulnerable  area  --  a plane  area  of 
suitable  size  and  shape  in  a plane  perpendicular  to  the  line  of 
sight  between  the  gun  system  and  the  missile  at  time  t; 

• a round  kills  the  cruise  missile  if  and  only  if  it  intercepts  the 
missile's  effective  vulnerable  area; 

• rounds  are  aimed  (at  time  t-tf)  at  the  estimated  center  of  the 
cruise  missile's  effective  vulnerable  area  at  time  t,  but  errors 
in  both  the  estimation  of  future  position  and  in  the  delivery 
points  of  rounds  may  occur. 

Under  these  conditions,  the  single-shot  kill  probability  SSKP(t)  at 
time  t may  be  modeled  as 

SSKP(t)  = A(t) *p(0,0)  (23) 

where: 

A(t)  = area  of  the  cruise  missile's  effective  vulnerable 
area  at  time  t,  and 
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p(0,0)  = value  of  a probability  density  function  of  the  two- 
dimensional  delivery  error  evaluated  at  the  aiming 
point. 

The  area  A(t)  of  the  effective  vulnerable  area  of  the  cruise  missile 
will  be  a function  of  the  dimensions  of  the  missile,  its  aspect  with 
respect  to  the  gun  system  at  time  t,  and  of  the  type  of  round  (fuze  or 
impact)  employed  by  the  gun.  In  detailed  simulations,  it  is  common  to 
model  targets  such  as  cruise  missiles  as  consisting  of  several  vulnerable 
components  with  vulnerable  areas  in  six  directions:  front,  left,  right, 
up,  down  and  rear,  relative  to  the  natural  coordinate  system  of  the  target. 
Areas  associated  with  each  direction  are  projected  first  onto  a plane 
perpendicular  to  the  closing  velocity  vector  between  projectile  and  target, 
and  then  onto  a plane  perpendicular  to  the  line  of  sight  between  weapon 
and  target.  It  has  been  found,  however,  that  simple  approximations  (such 
as  a spherical  representation)  to  the  multiple  component  model  produce 
results  which  are  often  satisfactory.  It  is  expected  that  the  spherical 
model  will  suffice  to  determine  the  effective  vulnerable  areas  A(t)  of 
cruise  missiles,  and  these  functions  are  therefore  regarded  as  being 
computable. 

Consider  next  the  intercept  probability  p(0,0).  For  each  round,  the 
delivery  error  at  the  time  of  predicted  intercept  may  be  modeled  as  a three 
dimensional  random  vector  <x2)in  which  each  component  x.  (i  = 1,2,3)  is 

_(x3) 

regarded  as  being  a sum£x-  . of  a number  of  random  variables  x-  - associ- 

• J * J 

ated  with  error  sources  which  contribute  to  errors  in  the  x.j -dimension. 
Under  the  assumptions  that 
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• for  each  i = 1,2,3,  the  x^.  are  independent, 

• for  any  pair  1 £ ij  f i2  1 i3»  the  x.  . and  x.  t are  independent 
for  all  j and  k,  and 

• the  number  of  error  sources  which  contribute  to  errors  in  each 
dimension  is  "large," 

(xi ) 

the  distribution  of  the  miss  vector  J x2  > may  (by  the  central  limit  theorem) 

Us) 

be  approximated  as  a trivariate  normal  distribution  without  correlation. 

To  utilize  the  effective  vulnerable  area  concept  described  above,  the  three 
dimensional  miss  vector  must  be  transformed  into  an  equivalent  vector  in 
the  plane  perpendicular  to  the  line  of  sight  between  target  and  weapon. 

This  transformation  may  be  accomplished  by  suppressing  the  range  dimension 
as  described  below.  As  will  be  seen,  this  transformation  produces  a bi- 
variate normal  distribution  (with  correlation)  for  the  equivalent  weapon 
delivery  error  vector  in  the  azimuth  elevation  plane. 

The  angle-like  errors  in  elevation  and  azimuth  are  obtained  by 
multiplying  angular  errors  by  range.  The  assumption  is  made  that  direct 
angular  errors  in  gun  pointing  lead  to  equivalent  angular  errors  in  the 
Cartesian  coordinate  system  perpendicular  to  the  line  of  sight.  Sources 
of  the  angular  errors  include  the  estimates  of  both  angles,  i.e.,  ele- 
vation and  azimuth,  and  estimates  of  the  first  derivatives  of  elevation 
and  azimuth  passed  from  the  tracking  system  to  the  fire  control  computer. 
These  errors  thus  lead  to  equivalent  errors  in  lead  angle  effects.  By 
assuming  that  errors  produced  by  the  fire  control  computer  are  negligible, 
and  by  neglecting  the  rotation  of  the  coordinate  system  with  range,  azi- 
much  and  elevation  axes  in  the  time  between  prediction  and  tracking,  it 
may  be  concluded  that  azimuth,  elevation  and  range  errors  produced  by  the 
tracking  system  produce  equal  errors  in  gun  pointing. 
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Range  errors  in  the  system  arise  from  two  primary  sources.  The  first 
of  these  is  the  tracking  system  which  provides  the  fire  control  computer 
with  estimates  of  range  and  the  first  derivative  of  range.  The  second 
source  of  range  error  is  the  estimation  of  muzzle  velocity.  This  error 
is  modeled  as  described  below. 

Suppose  that  the  projectile  at  time  of  predicted  intercept  has  veloc- 
ity W in  the  direction  (0,0,R)  and  that  the  miss  vector  at  time  of  predicted 
intercept  has  component  Z in  the  range  direction.  It  is  easily  seen  that 
the  time  required  for  the  projectile  to  reach  the  target  plane  is  approx- 
imately Z/( R-W)  where  R is  the  rate  of  change  of  range.  For  small  values 
of  R relative  to  the  value  of  W this  time  can  be  approximated  by  Z/W. 

With  respect  to  the  coordinate  system  with  axes  in  the  directions  of  range, 
elevation  and  azimuth,  and  by  neglecting  rotation  the  target  velocity  can 
be  expressed  as  (AE,  ER,  R)  where: 

A = azimuth; 

A = azimuth  rate  of  change; 

E = elevation; 

• 

E = elevation  rate  of  change; 

R = range; 

R = range  rate  of  change. 

The  error  component  Z thus  corresponds  to  errors  in  elevation  and  azimuth 
given  by  ERZ/W  and  ARZ/W  respectively. 

The  magnitude  of  the  range  error  is  the  sum  of  errors  due  to  the 
tracking  system,  which  are  assumed  to  pass  unaltered  through  the  fire  con- 
trol system,  and  errors  due  to  muzzle  velocity  variation.  The  latter  range 
error  is  modeled  by: 


28 


} 

i 

i 

I 


where: 


Range  error  = 


(Ji HI — (*\)  *. 

\1  + stf  (1  + stf)2  \d V/ / 


= projectile  flight  time; 

V = muzzle  velocity; 

s = slow-down  constant  of  projectile;  and 
dV  = muzzle  velocity  error. 


Under  the  assumptions  made  above,  probability  distributions  must  be 
obtained  for  the  following  errors: 

(1)  range; 

(2)  range  rate; 

(3)  azimuth; 

(4)  azimuth  rate; 

(5)  elevation  rate;  and 

(6)  muzzle  velocity. 

Since  normality  and  independence  have  been  assumed,  mean  or  bias  and  vari- 
ance are  sufficient  to  specify  the  distributions  required.  Note,  however, 
that  for  any  engagement  these  parameters  are  functions  of  time. 

The  following  form  for  the  probability  density  function  of  the  miss 
vector  in  a plane  perpendicular  to  the  line  of  sight  between  weapon  and 
target  incorporates  the  parameters  listed  above. 

Let  the  total  range-like  and  angle- like  errors  have  variances 


i 


given  by: 
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op  = variance  of  total  error  component  in  range; 

2 

ofl  = variance  of  total  error  component  in  azimuth;  and 
2 

oe  = variance  of  total  error  component  in  elevation. 
Then  transforming  the  range  error  into  the  azimuth-elevation  plane 
results  in  variances 

2 _ 2 . 2 A2  . 

°A  ' °a  +Jr  a"d 


2 . 2 . 2 E 

°E  °e  + °r  ^2  * 


where: 

2 

= variance  of  equivalent  azimuth  error; 

2 

ge  = variance  of  equivalent  elevation  error; 

A = azimuth  rate; 

E = elevation  rate;  and 
W = shell  velocity  at  time  of  intercept. 

Assuming  a two  dimensional  normal  distribution  and  neglecting  correlation, 
the  value  of  the  probability  density  function  at  the  target  center  is 


1 


27raEaA 


exp 


V 


B£2 

r~ 2 ' 7T 

2a£  2 


where : 


% 


B£  = bias  in  elevation; 

B^  = bias  in  azimuth  (i.e.,  the  mean  vector). 

One  final  factor  must  be  considered,  namely,  the  correction  necessary 
to  account  for  the  dependency  between  the  errors  associated  with  single 
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rounds  fired  in  a burst,  i.e.,  errors  common  to  all  rounds  in  the  burst. 
The  correction  is  obtained  through  consideration  of  the  burst  survival 
probability,  BSP,  given  the  burst  error  and  the  single  shot  survival 
probability  SSP.  Thus, 

BSP  = SSPn 

where  r.  is  number  of  rounds  in  the  burst. 

For  any  random  variable  X with  mean  x,  put  e * X - x.  Then 


E[eX]  = E[ex+e]  « E[ex+eex  + £ ex]  » 


- - .2  - 
2 


* ex  (1  + var(X))  . 


1 


Therefore, 


E[BSP]  = E[exp  (ntn  SSP)]  « exp  [ntnSSP]  (1  + - var  (tnSSP)). 

2 


Using  a Taylor  series  expansion  for  tn(l-x)  one  obtains 

2 

I 

~5~  + 

L\2o, 


— £n  SSP 


2iro^02 


exp 


•—  / n 2 

K]  "2 

— x—  + 


'1 


R 

2o 


)] 


where  R1  and  R2  are  the  components  of  the  burst  error  vector,  o*  and  a* 

are  the  intraburst  variances  and  A is  the  target  vulnerable  area.  Biases 

are  ignored.  Taking  expectations, 

A , 

-*nS$P  = • 1/2  1/2 

2ir(<j.j2+a32)  (a22+a^2) 

where  a32  and  of  are  the  variances  associated  with  the  bivariate  normal 
distribution  of  R^  and  R2.  A standard  calculation  yields 
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var  (tnSSP)  ■ (inSSP)2 


(°12^32)  (^22%2) 

ofz(2o3^)VHZ'W)'/2 


-1 


Implying  that  the  appropriate  correction  to  the  integrated  kill  rate  is 


This  correction  should  be  applied  over  each  burst  period;  It  should 
be  noted  that  if  the  higher  order  moments  of  tnSSP  are  large,  the 
approximations  used  here  are  no  longer  valid.  However,  this  is  rarely 
the  case  with  gun  systems. 

Assembling  the  methods  and  assumptions  outlined  above,  the  following 
formula  is  obtained  for  the  kill  rate  to  which  the  correction  factor 
for  inter-  and  intra-burst  errors  is  applied: 

FA(t)  / /B.  Bp\\ 

k(t)  = 7.1^17?  exp  1/2  W +hJ) 

where 

HA  = °A2R2  + °A2tf2R2  + 

+ *V  (»R2  + °R2  v)  • 

He  = oe2R2  + aE2tf2R2  + E2tf2R2 

+ K (aR2  + CTR2tf2 


5 
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and  where: 


JA1 


A2 

T 2 


0 2 

El 


E2 


r*2 


,*2 


R 


A(t) 

F 

• 

A 

E 

ctr2 

aR2 


= (baR+ca)2 
= (beR  + ce)2 

= °A12  + °A22 


2 - 


variance  azimuth  sensing 
variance  azimuth  gun-pointing 


“El2  + ° 


E2 


= variance  elevation  sensing 
= variance  elevation  gun -pointing 

= variance  azimuth  rate  input  to 
computer 

= variance  elevation  rate  input  to 
computer 

= range  (at  intercept) 

= time  of  flight 

= variance  in  muzzle  velocity 

= muzzle  velocity 
= vulnerable  area 

= average  firing  rate  during  a burst 
and  one  interburst  period 


(km2) 

(km2) 

(rad2) 

(rad2) 

(rad2) 

(rad2) 

(rad2) 

(rad2) 

/ rad2  \ 
\sec2  / 

/ rad2  \ 
\sec2  ) 

(km) 

(sec) 

(&) 

(fe) 

(km2) 
(sec-1 ) 


azimuth  rate  of  target  (at  intercept)  (rad/sec) 
elevation  rate  of  target  (at  intercept)  ^ 
variance  range  sensor  (input)  (km2) 

/ |(m2  \ 

variance  range  rate  input  to  computer  (sec7/ 
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ba 

= 

bal  + ba2  + ba3 

(rad) 

bal 

= 

sensor  azimuth  bias 

(rad) 

ba2 

= 

pointing  azimuth  bias 

(rad) 

ba3 

= 

A R"1  bR  tf 

(rad) 

b 

r 

= 

range  sensing  bias 

(km) 

be 

= 

bel  + be2  = be3 

(rad) 

bel 

= 

sensor  elevation  bias 

(radr 

be2 

= 

pointing  elevation  bias 

(rad) 

be3 

- 

;t  -R-^  :bR  t^ 

(rad) 

ca 

= 

cal  + CR  R 1 A bf 

(km) 

CR 

= 

prediction  error  due  to  nonlinear 
flight  in  range  direction1 

(km) 

cal 

= 

prediction  error  in  azimuth  direction1 

(km) 

ce 

= 

cel  + CR  R 1 E bf 

(km) 

cel 

= 

prediction  error  in  elevation 
direction1 

(km) 

For  the  most  part  the  parameters  in  the  above  list  must  be  supplied  by 
engineering  models  which  provide  values  for-  the  parameters  as  a function  of 
ship  systems  hardware  characteristics,  cruise  missile  characteristics,  and 
engagement  geometry  and  time.  DuHng  the  initial  stages,  of  the  project,  work 
was  carried  out  to  produce  engineering  models  which  predicted  biases  and 
variances  for  radar  systems.  This  work  is  described  in  appendix  D to 
[Bonder,  Cherry  and  Miller,  1973], 

'These  must  be  evaluated  as  a function  of  the  flight  path  in  any  specific  case 
under  study. 
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The  effects  of  EW,  in  particular  enemy  jamming  activities,  are  por- 
trayed in  the  model  through  the  effect*  of  these  activities  on  system  per- 
formance measures  in  the  above  list.  The  basis  of  engineering  models 
capable  of  predicting  the  relationships  between  certain  types  of  jamming 
and  performance  characteristics  of  radar  systems  is  described  in  appendices 
A and  B to  this  report. 

2.2.2  Kill  Rates  for  SAM  Systems 

Recall  that  the  kill  rate1  k(t)  for  a particular  SAM  system/cruise 
missile  pair  at  time  t was  defined  as 

6 - k(t ) = conditional  probability  that  the  cruise  missile  is 

killed  by  the  SAM  system  in  the  time  interval  [t,t+6), 
given  that  the  missile  is  still  alive  at  time  t and 
that  the  SAM  system  is  firing  at  the  cruise  missile 
at  that  time. 

The  probability  6 * k ( t ) that  the  cruise  missile  is  killed  by  a SAM  from  the 
SAM  system  in  the  interval  [t,t+6)  can  be  taken  in  the  form 

<5*K(t)  =• 1 (t)-SSKP  (24) 

where: 

i(t)  = probability  that  the  SAM  intercepts  the  cruise  missile 
in  the  time  interval  [t,t+s]; 


JAs  in  the  preceding  section,  the  notation'  «k(t)  is  used  instead  of 
•kij(t  d)  to  denote  a kill  rate  because  the  damage  state  and  the  weapon 
system  8 and  cruise  missile  j to  which  it  applies  is  understood. 


SSKP  = conditional  probability  that  the  SAM  kills  the  cruise 
missile,  given  that  it  intercepts  the  cruise  missile. 

As  the  above  notation  suggests,  and  in  contrast  to  the  situation  with 
rapid-fire  gun  systems,  the  single-shot  kill  probability  SSKP  for  SAMs  is 
here  assumed  to  be  a constant1.  The  value  of  this  constant,  which  may 
depend  on  the  type  of  cruise  missile  and.  the  type  of  SAM  but  not  on  time, 
may  be  estimated  using  existing  models.  It  is  regarded  as  being  a funda- 
mental descriptor  of  the  engagement.  The  determination  of  the  kill  rate 
6*k(t)  is  therefore  reduced  to  determining  i(t). 

Since,  by  assumption,  the  actual  flight  paths  and  speeds  of  the 
cruise  missile  with  respect  to  the  ship  are  deterministic,  the  flight  time 
the  SAM  requires  to  intercept  the  cruise  missile  at  any  point  along  its 
path  can  be  determined  for  the  particular  SAM  type  and  cruise  missile  type 
involved.  This  flight  time  then  determines  when  landing  must  have  taken 
place  in  order  for  the  SAM  to  intercept  the  cruise  missile  at  that  point 
on  its  flight  path.  One  may  therefore  write 

i(t)  = 6-J  h(t-u)  g(t,u)  du  (25) 


]This  assumption,  which  can  be  relatively  easily  relaxed,  seems  reasonable 
because  SAMs  are  (unlike  rounds  from  rapid-fire  guns)  guided  onto  targets. 
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where: 


h(x)»5  = Pr[SAM  fired  in  (x,  x + o)], 


and 

g(t,u)  = probability  density  function  for  the  flight  time  u 
of  a SAM  fired  to  intercept  the  cruise  missile  at 
time  t. 

Notice  that  the  above  development  of  i(t)  exhibits  the  dependence  of 
this  function  upon  the  past  in  that  the  probability  of  intercept  at  time 
t depends  upon  the  status  of  the  engagement  at  time  t - u.  Notice  also 
that  neither  of  the  functions  h(*)  and  g(t , - ) will  be  "smooth"  because 
the  firing  process  is  discrete  and  because  the  stochastic  variation  in 
flight  :.imes  for  the  SAM  will  normally  be  small.  However,  in  a limited 
sense,  this  model  removes  the  discrete  nature  of  the  arrivals  (or 
intercepts)  of  the  SAMs  at  the  cruise  missile. 

Under  the  assumption  that  the  SAM  flight  time  has  relatively  limited 
variation  --  say  a small -variance  uniform  random  variable  added  to  a 
deterministic  time  u(t)  — the  above  expression  for  i(t)  can  be  written 

u(t)+£ 

i(t)  = <5  • J h(t  - u)  2~  du 
u(t)-£ 


If  it  is  further  ..ssumed,  as  seems  reasonable,  that  the  probability 
6-h( t - u)  that  a SAM  is  fired  in  [t  - u,  t - u + 6]  has  the  constant 
value1 *  6*ht  for  u in  the  range  u(t)  - s,  <_  u(t)  + £,  the  above  expression 
becomes 


i(t)  = 6 • h. 


To  determine  the  values  of  the  ht  for  use  in  (26),  it  is  necessary 
to  consider  the  process  by  which  the  SAM  launching  system  operates.  The 
firing  process  from  the  SAM  launcher  can  be  described  by  a cyclic  Markov 
renewal  process  in  which  firing  takes  place  at  random  intervals.  As  in 
the  case  of  SAM  flight  times,  the  stochastic  variation  of  these  intervals 
may  be  small,  but  the  following  analysis  leads  to  a result  identical  to 
the  deterministic  case. 


1JS  the  notation  suggests,  the  "constant"  h.  may  have  different  values 

tor  different  times  t.  z 
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Suppose  that  there  are  R rails  on  the  SAM  launcher,  and  consider  an 
R state  Markov  renewal  process  with  transition  matrix: 
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where 

F-| (x)  = probability  distribution  of  time  to  reload. 

F(x)  = probability  distribution  of  time  between  firings  on  adjacent 
rails. 

Further,  let 
t-j  = mean  reload  time. 
t = mean  inter-firing  time. 

For  this  process,  it  can  be  shown  [Cinlar,  1969],  that  in  the  limit  the 
rate  at  which  firings  take  place  is  given  by: 

1 

T i + (R  - 1)t 


The  case  of  L independent  launchers  can  be  modeled  by  superposing  the 
models  for  the  individual  launchers.  In  this  case,  an  extension  of  re- 
sults in  [Cherry,  1972]  gives 

L 

t-j  + (R  - 1)t 

as  the  fate  at  which  firings  of  all  the  SAM  systems  take  place.  Using 
a technique  common  to  kill  rate  modeling,  the  probability  6*ht  that  a 
SAM  is  fired  in  the  interval  [t,  t + 6)  can  tnen  be  approximated  by  the 
mean  rate  at  which  SAMs  could  be  fired: 

Lt 

6 ' ht  ‘ ^ + (R  - 1)t  * 5 

where  Lt  is  the  number  of  launchers  available  at  time  t - u (u  being  the 
SAM  flight  time  yielding  intercept  at  time  t). 

Finally,  it  should  be  noted  that  the  above  discussion  has  dealt  with 

the  rate  at  which  SAMs  can  be  fired  whereas  the  equations  introduced  in 

the  preceding  section  deal  with  the  rate  at  which  SAM’s  arrive  at  the 
cruise  missile.  The  rate  at  which  SAMs  can  be  fired  does  not  correspond 
to  the  rate  at  which  intercepts  can  occur  since  a Doppler  effect  is 
present  due  to  the  flight  times  of  the  SAMs  and  the  motion  of  the  tar- 
get cruise  missile.  Let 

t-j  = time  of  first  intercept, 

t^  = time  of  second  intercept, 

T.|  = time  of  flight  of  first  SAM,  and 

T2  = time  of  flight  of  second  SAM. 
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Then  the  time  between  launchings  is  given  by; 

<‘2  - T2>  - (t,  - T,) 

* (t2  - t1)  - (T2  - T1) 

The  length  of  the  inter-intercept  time  period  in  the  limit  as  the  period 
grows  small  is  equal  to  the  firing  interval  times  1 - DtT^  where  T^  is 
the  time  of  flight  and  Dt  indicates  differentiation  with  respect  to  in- 
tercept time.  The  model  firing  rate  thus  becomes: 

0 - DtTf) 

Note  that  the  above  rate  corresponds  to  the  situation  in  which 
guidance  channel  constraints  are  not  operative.  In  the  case  in  which 
guidance  channel  constraints  are  operative  the  rate  at  which  SAMs  can 
be  fired  corresponds  to  the  rate  at  which  intercepts  occur.  Consider  the 
non-homogeneous  Poisson  process  with: 

- J k(u)du 

Pr[T  > t]  = e 0 

For  this  process  it  can  be  shown  that  the  mean  rate  of  events  over  a time 
period  [O.t]  is  given  by 


I 


Accordingly,  approximate  the  mean  rate  of  kills  by  the  expression: 
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M 

f,(t>  - -g 

j=i 


in  m.(t) 
J- 


where 


f,(t) 

mj(t) 


M = 


mean  rate  of  kills  in  the  period  [0,t], 
Pr[cruise  missile  j survives  to  time  t], 
number  of  cruise  missiles  engaged. 


This  expression  reflects  the  number  of  kills  per  time  unit;  the  rate 
required  is  the  number  of  intercepts.  Consider  a Poisson  process  with 
parameter  x and  suppose  that  each  time  an  event  occurs  in  this  process 
it  is  recorded  with  probability  SSKP.  It  can  be  shown  that  the  recorded 
event  process  is  Poisson  with  parameter  X*SSKP.  Based  on  this  analogy, 
the  approximation  f(t)  = SSKP"^  • f-j(t)  is  used  for  the  firing  rate  under 
guidance  constraints.  The  firing  rate  used  thus  becomes 

f(t)  (1  - DtTf) 

for  those  cases  in  which  f(t)  is  less  than  the  number  of  guidance  channels 
available. 


2.3  Determination  of  Firing  Probabilities 

The  firing  probabilities  p..(t|d),  which  measure  the  likelihood  that 

' J 

each  weapon  system  i is  actually  firing  at  each  cruise  missile  j at  time 
t,  may  be  modeled  as  a function  of: 

• whether  or  not  cruise  missile  j has  already  been  destroyed 
or  has  impacted  the  ship  by  time  t;  and 

• the  logic  whereby  the  ship  assigns  defensive  weapons  to 


surviving  cruise  missiles. 
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It  has  been  assumed  (see  the  list  of  assumptions  near  the  beginning  of 
section  2.1)  that  the  ship  uses  a weapon  assignment  logic  in  which  all 
available  defensive  weapons  are  always  assigned  to  the  one  cruise  missile 
which  will,  if  not  sooner  destroyed,  be  the  next  to  impact  the  ship.  This 
logic  requires  that  the  ship  has  the  capability  to  evaluate  the  relative 
threats  posed  by  the  surviving  cruise  missiles  at  any  time.  A conceptual 
threat  evaluation  system  that  the  ship  might  use  for  this  purpose  is 
therefore  briefly  discussed  in  section  2.3.1.  The  manner  in  which  this 
threat  evaluation  system  could  provide  data  for  the  assumed  weapon 
assignment  logic  is  discussed  in  section  2.3.2.  Finally,  expressions 
for  the  firing  probabilities  which  reflect  the  structure  of  the  assumed 
weapon  assignment  logic  are  developed  in  section  2.3.3. 


2.3.1  Threat  Evaluation 


For  a single  ship/multiple  cruise  missile  engagement,  the  relative 
threat  posed  by  each  of  several  incoming  cruise  missiles  may  be 
evaluated  in  terms  of:1 

(1)  the  predicted  length  of  time  until  impact;2 

(2)  the  extent  to  which  defensive  weapon  systems  are  already 
committed  to  handling  previously  evaluated  cruise  missiles;  and 


xIn  the  case  of  a multiple  ship/multiple  cruise  missile  engagement, 
other  factors,  such  as  the  value  of  each  ship,  should  be  added  to  the 
list  of  items  in  terms  of  which  threats  should  be  evaluated. 

depending  on  the  configuration  and  capabilities  of  the  defensive  weapons 
on  the  ship,  the  position  of  a cruise  missile  at  each  time  may  also 
be  important  to  an  evaluation  of  the  threat  posed  by  the  missile  at  that 
time.  For  example,  a cruise  missile  may  be  observed  to  be  about  to  enter 
a portion  of  its  track  which  is  not  coverable  by  the  ship's  defensive 
weapons  because  of,  say,  elevation  or  azimuth  limitations  of  these  wea- 
pons so  that  the  cruise  missile,  if  not  immediately  engaged,  will  surely 
impact  the  ship. 


(3)  the  extent  to  which  defensive  actions  have  been  taken  whose 
outcomes  have  not  yet  occurred. 

The  length  of  time  remaining  until  impact  is  clearly  relevant  to 
evaluating  the  threat  posed  by  an  incoming  cruise  missile.  Were  it  not 
for  items  (2)  and  (3)  in  the  above  list,  the  threats  posed  by  incoming 
cruise  missiles  could  perhaps  reasonably  be  ranked  solely  on  the  basis  of 
remaining  time  to  impact.  Indeed,  even  when  factors  (2)  and  (3)  are 
acknowledged,  one  may  still  argue  that  ranking  threats  by  remaining  times 
to  impact  regardless  of  the  defensive  actions  which  may  have  been  taken 
against  some  of  the  missiles  and  whose  outcomes  have  not  yet  occurred 
may  be  a conservative  strategy.  As  (2)  and  (3)  suggest,  however,  some 
schemes  for  threat  evaluation  are  structured  to  account  for  the  fact  that 
resources  already  committed  to  a cruise  missile  may,  after  a time  lag, 
result  in  the  destruction  of  that  missile  so  that  no  further  action  against 
the  missile  would  be  necessary.1  However,  the  previously-made  assumption 
that  the  ship  can  redirect  in-flight  fire  as  long  as  the  fire  has  not 
reached  the  cruise  missile  at  which  it  was  last  targeted  effectively  as- 
sumes away  items  (2)  and  (3)  in  the  above  list.  It  is  therefore  consis- 
tent to  assume  that  the  ship  uses  a threat  evaluation  system  in  which  the 
threat-rank  of  each  cruise  missile  is  the  same  as  the  order  in  which  it 
will,  if  not  sooner  destroyed,  impact  the  ship  and  in  which  the  threat 


]Threat  evaluation  structures  of  this  second  type  could  be  designed  by 
ranking  threats  by  weighted  remaining  times  to  impact.  The  weights 
assigned  to  cruise  missiles  would  be  functions  of  the  probabilities 
that  these  missiles  will  survive  the  defensive  actions  already  taken. 
Threat  evaluations  would  then  be  accomplished  in  a computational  envi- 
ronment involving,  in  addition  to  the  simple  flight  time  calculations, 
calculations  of  a more  complex  character. 
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rank  of  a cruise  missile,  once  assigned,  changes  only  when  the  missile  is 
destroyed  or  impacts  the  ship1  and  does  not  otherwise  change  in  time. 

This  assumed  threat  ranking  procedure  will  be  used  to  provide  input  to 
the  assumed  ship's  weapon  assignment  logic  which  is  discussed  in  the  next 
section. 

2.3.2  Weapon  Assignment 

In  a single  ship/multiple  cruise  missile  engagement,  the  objective 
of  a weapon  assignment  logic  is  to  assign  defensive  weapons  to  incoming 
cruise  missiles  in  such  a way  as  to  minimize  some  measure  of  the  damage 
to  the  ship  (e.g.,  expected  "damage",  probability  of  impact  of  one  or 
more  cruise  missiles,  etc.).  To  accomplish  this  objective,  a weapon 
assignment  procedure  should  account,  for  at  least  the  following: 

(1)  the  relative  threat  posed  by  each  of  several  incoming 
cruise  missiles; 

(2)  the  availability  of  defensive  weapon  systems;  and 

(3)  the  characteristics  of  the  defensive  weapon  systems  and 
of  the  cruise  missiles. 

The  relative  threats  posed  by  the  incoming  cruise  missiles  can  be 
determined  by  a suitable  threat  evaluation  scheme  as  discussed  previously. 
The  notion  of  defensive  weapon  system  availability , as  used  above,  in- 
cludes a consideration  of  whether  or  not  a given  defensive  weapon  (to- 
gether with  its  required  support  systems  such  as  guidance  links,  gun 


bruise  missiles  which  have  been  destroyed  or  which  have  impacted  the 
ship  may  be  said  to  have  a null  threat  rank. 


t 
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pointing  radars,  etc.)  is  still  alive  as  well  as  a consideration  of  whether 
or  not  the  weapon  is  already  in  use  and,  if  so,  whether  or  not  reassign- 
ment might  be  profitable.  The  characteristics  of  the  defensive  weapon 
systems  which  are  relevant  to  the  weapon  assignment  process  clearly  in- 
clude such  factors  as  range,  accuracy,  rate  of  fire,  time  of  flight  for 
projectiles,  etc. 

Two  weapon  assignment  logics  have  been  described  by  [Forsyth  etal., 
1973]1.  The  first  of  these  procedures  has  the  advantage  of  being  ex- 
tremely simple.  However,  it  appears  that  this  procedure  may,  when  em- 
ployed against  cruise  missiles,  result  in  random  weapon  assignments.  The 
second  procedure  is  considerably  more  complex;  it  involves:2 

•the  ranking  of  incoming  cruise  missiles  by  threat  (the 
threat  ranks  may  change  in  time  and  may  depend  on  defensive 
actions  already  taken); 

•the  establishment  of  categories  of  threat  ranks; 

•the  assignment  of  defensive  weapons  to  cruise  missiles 
whose  threat  rank  is  in  the  highest  threat  category  until 
this  category  becomes  empty  (this  may  involve  the 
reassignment  of  defensive  weapons  assigned  to  cruise  missiles 
with  threats  in  a lower  category  to  cruise  missiles  whose 
threats  enter  a previously  empty  higher  category). 


^oth  of  the  weapon  assignment  procedures  described  by  Forsyth  include 
some  aspects  of  what  has  here  been  called  "threat  evaluation"  and  been 
viewed  as  being  distinct  from  weapon  assignment. 

2The  structure  of  this  second  weapon  assignment  procedure  resembles  that 
of  a "foreground/background  queue",  and  advantage  of  this  correspondence 
could  be  taken  in  constructing  a reasonably  detailed  model  of  the  pro- 
cedure. 


TV 


The  second  of  the  above  weapon  assignment  procedures  appears  to  be  the 
more  promising  of  the  two. 

In  view  of  assumptions  which  have  previously  been  made,  namely, 

• that  the  ship  uses  a threat  evaluation  system  in  which 
threat  ranks  of  a cruise  missile,  once  assigned,  change 
only  when  the  missile  is  destroyed  or  impacts  the  ship 
(see  the  discussion  in  the  previous  section),  and 

• that  the  ship  has  the  capability  to  redirect  in-flight 
SAMs  and/or  gunfire  to  new  targets  as  long  as  these 
projectiles  have  not  reached  the  cruise  missile  at  which 
they  were  last  targeted  (see  the  discussion  at  the  beginning 
of  section  2,0), 

the  assumed  weapon  assignment  logic  is  a special  case  of  the  second  weapon 
assignment  procedure  described  by  Forsyth  in  which  the  threat  categories 
of  the  latter  consist  of  singleton  sets. 

2.3.3  Determination  of  Firing  Probabilities 

The  assumptions  which  have  been  made  in  the  preceding  two  subsections 
about  the  ship's  threat  evaluation  system  and  weapon  assignment  logic  may 
be  summarized  as  follows: 

• the  threat  rank  assigned  by  the  ship  to  a cruise  missile  is  the 
same  as  the  order  in  which  the  cruise  missiles  will,  if  not  sooner 
destroyed,  impact  the  ship; 

• the  threat  rank  of  a cruise  missile,  once  assigned,  changes  only 
when  the  missile  is  destroyed  or  impacts  the  ship  and  does  not 
otherwise  change  in  time; 
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• at  any  time  during  the  engagement,  all  available  defensive 
weapons  are  assigned  to  the  cruise  missile  having  the  highest 
threat  rank  and  remain  so  assigned  until  the  missile  is  destroyed 
or  impacts  the  ship. 

Recall  that,  by  definition,  p..(t|d)  is  the  conditional  probability 

• J 

that  defensive  weapon  system  i is  assigned  to  (and  firing  at)  cruise  missile 
j at  time  t,  given  that  the  ship  is  in  damage  state  d and  that  cruise 

missile  j is  still  alive  at  that  time.  In  view  of  the  assumed  threat 

evaluation/weapon  assignment  logic  described  above,  either  all  the  ship's 
available  defensive  weapons  will  be  assigned  to  a particular  cruise  missile 
at  any  time,  or  none  will  be  so  assigned.  In  particular,  the  available 
defensive  weapons  will  all  be  assigned  to  cruise  missile  j at  time  t if 
and  only  if 

• cruise  missile  j is  still  alive  at  time  t,  and 

• either  j=l  or,  if  j >2,  cruise  missile  1,  ...,  j-1  have 

beer,  destroyed  or  have  impacted  the  ship  by  time  t. 

Thus 


n>1  (t  | u) 


P "j  ( 1 1 d ) = 


undefined 

0 


if  d=u  and  0 < t < t-j , 
if  d^u  and  0 <_  t < t^ , 
otherwise. 


(27) 


Moreover,  we  have  on  the  one  hand  p..(t|d)  = 0 for  all  1 < i < W,  all  d, 

* J 

all  2 £ j <_  M and  all  t > t..  On  the  other  hand,  for  all  1 <_  i <_  W,  all  d, 

J 

all  2 <_  j <_  M and  all  t < t - , we  have 

J 


p . . ( 1 1 d)  = Pr[i  firing  on  j at  time  tjship  in  d and  j alive  at 

' J 

time  t] 
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Pr[i  firing  on  j arid  j alive  at  t[ ship  in  d at  t] 
Pr[j  alive  at  t | ship  in  d at  tj 


PrCl » •••,  j-1  dfeud  at  t and  j alive  at  t|ship  d at  t] 
Pr[j  alive  at  t | ship  in  d at 


Pr[j-1  dead  at  t and  j alive  at  t | ship  in  d at  t] 
Pr[j  alive  at  t|ship  in  d at  t] 


1 

Pr[j  alive  at  t J ship  in  d at  t] 


»|jPr[j-l  dead  at  t and  j alive  at  t | ship  in  d at  t] 

+Pr[j~l  alive  at  t and  j alive  at  t J shi p in  d at  t]} 

-Pr[j-1  alive  at  t and  j alive  at  t | shi p in  d at  t]J 


1 


Pr[,j  alive  at  tjship  in  d at  t] 


|pr[j  alive  at  t | ship  in  d at  t] 

-Pr[j-1  alive  at  tjship  in  d at  t]j 


1 


sritw  ‘ c”j(tid)  ■ Vi (t|d)] 

J 


= [1  -m._1(t|d)/m.(t|d)]. 


(28) 


The  third  equality  in  the  above  string  of  equations  follows  from  the  fact 
that,  according  to  the  weapon  assignment  logic,  defensive  weapon  system  i 
will  be  firing  at  cruise  missile  j >_2  at  time  t if  and  only  if  cruise 
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missiles  1,  •••,  j-1  are  dead  at  that  time  and  cruise  missile  j is  alive 
at  that  time.  The  fourth  equality  follows  from  the  fact  that  the  death 
of  cruise  missile  j-1  implies  that  of  the  preceding  cruise  missiles,  if 
any.  The  sixth  equality  follows  from  the  fact  that  the  aliveness  of  cruise 
missile  j-1  implies  that  of  cruise  missile  j (because  of  the  structure  of 
the  assumed  weapon  assignment  logic).  The  other  equations  are  merely 
manipulative;  they  express  facts  about  conditional  probabilities  or,  in 
the  case  of  the  last  two  equations,  the  definition  of  the  m- ( • | • ) . 

Equations  (27)  and  (28)  may  be  used  in  (21)  to  compute  the  survival 
probabilities  m.(*|d)  for  the  cruise  missiles  and  the  damage  state 

J 

probabilities  qd(*)  for  the  ship  for  the  case  when  the  ship  uses  the  threat 

evaluation/weapon  assignment  logic  which  has  been  assumed.  As  indicated  in 
in  the  discussion  which  led  to  (21),  the  normal  procedure  would  be  to 

bootstrap  (21)  numerically  to  obtain  the  survival  probabilities  m-(t|d) 

J 

and  p..(t|d)  for  all  time  t in  each  time  interval  [t^,  tk+^).  However, 
the  simple  form  of  (27)  and  (28)  suggests  that  it  may  be  possible  to 
insert  these  expressions  in  (21)  and  integrate  the  resulting  system  of 
equations  directly,  thus  bypassing  the  numerical  bootstrapping.  This 
approach  has  been  examined.  While  it  appears  that  the  use  of  numerical 
techniques  cannot  be  completely  avoided  by  using  the  analytic  approach, 
it  is  possible  to  reduce  the  numerical  problem  to  a form  to  which  simple 
and  well-known  numerical  methods  can  be  applied.  Details  are  given  in 
appendix  C. 


Preceding  page  blank 


3.0  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


The  single  ship/multiple  cruise  missile  engagement  model  described 
in  the  previous  section  has  been  developed  to  a point  where  computer  im- 
plementation and  initial  tryout  is  now  appropriate.  VRI  feels  that  the 
model  has  considerable  potential  as  an  aid  to  analysis  of  ship/cruise 
missile  engagements  and  that  the  model  may  be  helpful  in  devising  counter- 
measures (e.g.,  improved  or  satisfactory  weapon  lethality  characteristics 
with  respect  to  cruise  missiles)  to  the  cruise  missile  threat.  VRI  there- 
fore recommends  that  computer  implementation  and  initial  tryout  be  scheduled. 

VRI  also  recommends  that  the  initial  runs  of  the  model  include  a varia- 
tional analysis  to  establish  the  sensitivity  of  the  model  outputs  to  the 
inputs.  A high  degree  of  sensitivity  of  the  outputs  to  the  inputs  will 
suggest  the  need  for  accuracy  in  the  input  data  so  as  to  ensure  reliable 
outputs. 

In  addition  to  establishing  the  extent  to  which  submodeling  may  be 
necessary  for  data  generation,  the  initial  trial  runs  of  the  model  may 
suggest  that  modifications  to  the  existing  model  structure  are  appropriate. 
While  it  is  unlikely  that  changes  to  the  overall  mathematical  structure  of 
the  model  would  be  indicated  as  a result  of  the  trial  runs,  it  must  never- 
theless be  recognized  that  the  development  of  most  useful  models  is  an 
evolutionary  process  and  that  structural  changes  made  in  the  light  of  ex- 
perience gained  in  using  the  model  cannot  be  ruled  out. 

Aside  from  the  research  which  may  be  necessary  for  data  generation  or 
evolutionary  structural  modification  as  mentioned  above  (the  need  for  either 
of  which  is  not  yet  established),  there  appear  to  be  two  types  of  research 


which  may  be  necessary  in  order  for  the  model  to  realize  its  full  potential; 
these  are: 

• refinements  - improvement  in  the  flexibility  or  degree  of 

realism  with  which  the  model  portrays  engage- 
ments in  the  scenario  as  presently  conceived;  and 

♦ extensions  - enlargements  of  the  scenario  that  the  model 

portrays. 

The  most  obvious  example  of  a refinement  to  the  existing  model  --  and 
one  which  should  be  a high-priority  objective  of  future  research  — is  the 
portrayal  of  different  weapon  assignment  logics  for  the  ship.  Indeed,  the 
ship's  weapon  assignment  logic  is  one  of  the  primary  variables  of  the  en- 
gagement which  is  subject  to  the  control  of  the  ship,  and  so  the  degree  to 
which  the  model  (as  presently  structured  or  as  modified  in  the  future)  will 
be  useful  in  devising  countermeasures  to  the  cruise  missile  threat  will 
therefore  be  in  direct  proportion  to  the  degree  of  flexibility  the  model 
has  in  portraying  different  and  varied  weapon  assignment  logics. 

A second  refinement  to  the  existing  model  is  the  relaxation  of  the 
assumption  that  the  ship  can  redirect  in-flight  gunfire  and  SAMs;  in 
VRI's  opinion,  however,  this  is  definitely  secondary  in  importance  when 
compared  to  endowing  the  model  with  the  capability  to  portray  different 
and  varied  weapon  assignment  logics. 

The  most  obvious  example  of  an  extension  to  the  existing  single  ship/ 
multiple  cruise  missile  engagement  model  --  and  one  that  should  receive 
serious  attention  in  the  near  term  --  is  the  enlargement  of  the  existing 
model's  scenario  to  include  more  than  one  ship.  The  result  of  research  of 
this  type  would  be  a "multiple  ship/multiple  cruise  missile  engagement 
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model"  which  could  be  useful  as  an  aid  to  planning  for  fleet  air  defense 
against  cruise  missiles.  It  may  be  expected  that  such  a multiple-ship 
model  will  at  least  superficially  resemble  the  single-ship  model  and  may 
in  fact  have  a closely  parallel  structure,  so  that  the  single-snip  model 
may  be  regarded  as  being  a prototype  of  the  multiple-ship  model  — a 
natural  step  in  an  evolutionary  chain  of  models.  The  fact  that  present 
Navy  doctrine  often  calls  for  ships  (some  of  enormous  value  in  terms  of 
construction  cost)  to  travel  in  company,  together  with  the  serious  and 
recognized  threat  posed  by  cruise  missiles,  suggests  that  the  multiple 
ship  model  be  developed  as  soon  as  practicable. 
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APPENDIX  A 

RADAR  ANGULAR  AND  RANGE  TRACKING  ERRORS: 

AMPLITUDE  COMPARISON,  SEQUENTIAL  LOSING  AND  ANGLE  TRACKING 

A submodel  for  predicting  values  of  the  angular  tracking  errors  for 

2 2 

gun-system  radars  (these  are  the  quantities  denoted  by  o^l  and  o^2  in 
section  2.2. i)  in  the  presence  of  enemy  jamming  is  described  in  this 
appendix.  Theoretical  lower  bounds  on  the  variances  of  angle  error 
estimators  are  derived,  and  some  of  the  most  commonly  used  error  estima- 
tion techniques  are  analyzed  to  determine  how  closely  they  approach 

A 

these  theoretical  lower  limits  (Cramer  Rao  bounds). 

A 

More  specifically,  this  appendix  derives  the  Cramer  Rao  bound  on  the 
variance  of  error  angle  estimators  and  then  discusses  in  some  detail  one 
implementation  of  such  an  estimator  called  amplitude  comparison,  sequen- 
tial lobing.  The  first  and  second  moments  of  the  conventional  sequential 

A 

lober  are  derived  and  comparisons  are  then  made  with  the  Cramer- Rao  bound. 
The  derivation  is  largely  based  on  an  unpublished  report1  issued  by  Tech- 
nology Service  Corporation  [Lank,  Poll  on,  1969]. 

Appendix  B will  discuss  a second  technique  called  amplitude  compari- 
son monopulse,  which  is  also  used  to  obtain  angle  error  information. 

The  results  concerning  angle  tracking  errors  can  be  directly  used 
as  an  input  to  a target  tracking  analysis.  See,  for  example,  appendix 
D in  [Bonder,  Cherry,  Miller,  1973]. 

Complex  signal  notation,  as  shown  for  example  in  [Miller,  1969]  and 
[Reed,  1962],  is  employed.  The  notation  E(*)  stands  for  the  expectation 


^he  report  appears  in  IEEE  Trans.  Aerospace  and  Elect.  Sys.  , Vol.  AES-10, 
No.  3,  May  1974,  pp.  393-397. 


of  the  quantity  in  the  parenthesis  while  Re(*)  and  Im( • ) stand  for  the 
real  and  imaginary  parts  of  the  indicated  argument.  A superscript  * 
denotes  the  conjugate  of  a complex  number,  and  absolute  values  are 
denoted  by  |*|.  A superscript  p denotes  the  transpose  conjugate  of  a 
vector  or  matrix.1 

A.l  Cramer-Rao  Bound:  Noncoherent  Processing 

It  is  assumed  that  M noncoherent  received  waveforms  or  pulses  are 
to  be  processed  in  order  to  obtain  the  direction  and  the  magnitude  of  the 
error  angle,  denoted  by  e,  between  the  target  direction  and  the  zero  track- 
ing error  direction  (which  is  usually  the  axis  of  symmetry  of  the  antenna) 
in  one  dimension  or  coordinate.  The  purpose  of  this  section  is  to  derive 
the  Cramer-Rao  lower  bound  on  the  variance  of  an  unbiased  estimator  of 
which  utilizes  the  following  complex  observables  when  noncoherent  pro- 
cessing is  employed: 


Zm  = AGm(e)e  m + V m = (A1> 


^ther  notation  will  be  introduced  as  needed.  All  notation  used  in  this 
appendix  is  independent  of  that  used  elsewhere  in  this  report  even  though 
some  overlap  occurs. 
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where 

e is  the  error  angle  between  the  target  direction  and 

the  zero  tracking  error  direction, 

/ , 1 

G(-)  are  the  normalized  two-way  antenna  yol tags  patterns 

measured  relative  to  the  zero  tracking  error  direction, 

A is  the  observable  amplitude  which  is  normalized  with 

respect  to  the  antenna  beam  boresight  direction. 

th  1 

G is  the  phase  of  the  m observable, 
m 

M is  the  number  of  noncoherent  received  waveforms  that  are 
processed. 

The  are  t'ne  signal's  caused  by  system  noise  and  are  zero  mean,  complex 

Gaussian  random  variables  which  are  assumed  mutually  independent.  In 

addition,  the  in-phase  and  quadrature  components  (real  and  complex  parts) 

2 

of  each  are  assumed  independent  with  the  same  variance  u . Under  these 
conditions,  the  complex  Gaussian  random  variables  have  the  properties 
[Reed,  1962] 


E<v„0 


’The  receive  pattern  is  not  necessarily  the  same  as  the  transmit  pattern 
and  distinct  observables  may  be  received  with  the  same  antenna  pattern. 

2The  case  when  the  phase  is  the  same  for  each  target  signal,  i.e.,  monopulse 
processing,  will  be  discussed  in  appendix  B. 
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and  all  other  first,  second,  third,  and  fourth  order  moments  are  zero,  where 
2 

a is  the  variance  of  the  in-phase  and  quadrature  noise  variates, 

* is  the  complex  conjugate  operation. 

The  complex  probability  density  function  of  the  Z^  can  be  written  as 

p = = C exp  [-  \ (Z  - TOY1  (Z  - H)1  , (A2) 

where 


C is  a normalizing  constant, 

K = j E(VVP)  = o2I  , 

I is  the  identity  matrix, 

p is  the  conjugate-transpose  operation. 

There  are  M + 2 unknown  parameters:  e,  A,  e^,  m = !,..., M. 


63 


The  likelihood  function  can  now  be  written  as 


L = an  p = An  C - — ZPZ  - 2Re(ZpH)  + A2  TV  (c) 

L — ' n _ ' 


wnere 


Re  is  the  real  part, 

all  summations  go  from  1 through  M. 

Expanding  the  middle  term  of  (A3)  it  is  seen  that  the  likelihood  function 
can  be  written  as 


L = £n 


C \ ZPZ  - 2A  yT*  G (e)(x  cos  6^  + y sin  8 ) + A2  G2(e) 

2 2 [_  / j m ' m m m L — / nr  _J 


Making  use  of  (A4)  it  is  seen  that  the  (M  + 2)  x (M  + 2)  Fisher  in- 
formation matrix  has  the  form 


E (i  ui 


-E  ( -Mr 


04 


a12  = a21 


3 l\  a *,  % 

s^r) = s?  Zj 


l,k  + 2 


2,1  = • °»  J = 


2,k  + 2 


'Efe)  = °* 


aj  + 2,2  = °>  J = ^ 


Thus  the  information  matrix  N can  be  written 


P 0 
0 Q 


a11  a12 


a21  a22 
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Q is  a H x M square  matrix, 
all  other  elements  of  N are  zero. 

If  e is  an  unbiased  estimator  of  e,  the  error  angle  between  the  targe 

2 

direction  and  the  zero  tracking  error  direction,  with  variance  o , then 

e 

2 

the  Cramer-Rao  theorem  states  that  a lower  bound  for  a , when  noncoherent 

e 

processing  is  utilized,  has  the  form 


U, 


(A6) 


where  Ui  is  a vector  having  a one  as  its  first  element  and  zeros 
elsewhere.  Because  of  the  form  of  N given  by  (A5)  this  can  be  written  as 


U 1 P \ 


(A7) 


Making  use  of  the  expressions  developed  above  it  is  seen  that  when 
noncoherent  received  waveforms  are  processed  (A7)  can  be  written  as 


E M 


- 2R 


_E  (GmU>)  __E  _ ' _EGB.U)Gra(e)_ 


(AS) 


rr 


where  R is  the  per  pulse  signal-to-noise  ratio  for  each  beam  and  it  is 
assumed  that  the  target  is  located  at  the  boresight  of  the  beam 


R 


> 


Es  is  the  energy  in  each  received  waveform  and  Nq  is  the  one-sided  power 
spectral  density  of  the  system  noise. 

For  M = 2,  which  is  the  case  when  two  antenna  patterns  are  used  to 
derive  the  error  angle  e,  it  is  seen  that  (A8)  can  be  written  as. 


1 

2R(G*  + g|) 


(A9) 


where 


= G-j(e)  , = G^{e) 


Equation  (A9)  is  the  one  dimensional  analog  of  equation  (69),  page  29, 
in  [Horstetter&  Delong,  1969]. 
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A.?.  Cramer-Bao  Bound:  Coheri.  Processing 


If  we  use  the  same  notation  as  developed  in  the  previous  section, 
then  the  complex  radar  observables  when  coherent  processing  is  employed 
have  the  form 


Z = AG  (eje1^-0)  + V , m=l,...,M 

in  m * n 


where 


is  the  doppler  frequency  (assumed  very  much  smaller  than  the 
transmitted  center  frequency), 

are  known  time  delays  relating  the  times  of  occurrence  of  the  Z^. 

For  convenience  it  is  assumed  that  = 0. 

is  the  unknown  phase  corresponding  to  the  first  observable. 


Thus  there  are  four  unknown  parameters:  e.  A,  w,  e. 

The  likelihood  function  corresponding  to  (A10)  can  now  be  written  as 


L = c ' [z°z  - 2 cos  (UtD  - 6)  + yB  sin  - e)j 


The  elements  of  the  4 x 4 Fisher  information  matrix 


" ■ (aij> 
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have  the  following  values 


a 


11 


» 


5 


A 

2 


]Cg®(e1g'»(':) 


al3  = a31  = 


= 0 


aK  = a41 


y 


-E 


0 
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? 

a23  = a32 


a24  = a42 

Thus  by  using  the  same  procedure  as  employed  in  the  previous  section  it  can 
be  shown  that  the  lower  bound  of  the  variance  of  an  unbiased  estimator 
of  e,  the  error  angle  between  the  target  direction  and  the  zero  tracking 

error  direction,  is  also  given  by  ( A8  ) for  the  case  of  coherent  pro- 
cessing. Therefore,  the  Cramer-Rao  lower  bound  is  the  same  if  either  non- 
coherent and  coherent  processing  is  used  to  obtain  an  estimate  of  the  error  angle. 

If  the  signal-to-noise  ratio  is  high  enough,  then  well-designed 
processors  should  achieve  about  the  same  performance,  in  terms  of  estimator 
variance,  for  either  non-coherent  or  coherent  radar  systems.  This  point  is 
noted,  for  example,  on  page  28  of  the  paper  by  [Hofstetter  and  DeLong,  1969], 
Therefore,  from  this  point  on  we  shall  only  discuss  angle  estimation  for 
noncoherent  radar  systems. 


/u 


A..->  I)ir./‘nr.r.f.on  of  t'.u>  Tion-Br.im  lh'.t ! motor 

Expression  (A9)  gives  the  Cramer-Rao  lov/er  bound  on  the  variance 
of  an  unbiased  estimator  of  the  tracking  error  angle  e when  two  received 
waveforms  are  used  to  derive  the  tracking  error  signal.  If  M received 
waveforms  are  processed  for  each  of  two  antenna  beam  positions  then  it  is 
easily  shown  that  the  expression  for  the  lower  bound  must  be  modified  to 
read 


2 

o 

e 


(A12) 


where  the  notation  given  in  the  previous  sections  is  used.  Almost  all 
operational  tracking  error  estimators  use  the  received  signals  from 
two  antenna  beam  positions  to  derive  tracking  error  information  and  the 
remainder  of  this  paper  will  only  be  concerned  with  such  two  beam  estimators. 
From  this  point  on  it  will  also  be  assumed  that  each  antenna  beam  is 
symmetric  about  its  boresight  direction  and  that  the  beams  are  squinted 
or  offset  by  the  same  amount  on  either  side  of  the  zero  tracking  error 
direction  (which  is  usually  the  symmetry  axis  of  the  antenna).  The  same 
antenna  pattern  (referenced  to  the  boresight  direction)  is  assumed  for 
each  of  the  beams.  The  angle  between  the  zero-tracking  error  direction 
and  the  antenna  beam  boresight  directions  will  be  denoted  by  6.  A picture 
of  the  geometry  involved  is  given  in  figure  Al. 
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• All  angles  are  measured  in  the  counter-clockwise  direction 

• p is  the  tracking  error  angle 

• <5  is  the  antenna  beam  squint  or  offset  angle 

• ?(•)  is  the  normalized  two-way  antenna  voltage  pattern  measured 
with  respect  to  the  beam  boresigbt 

• G-|(*)  and  G^i’)  are  the  normalized  two-way  antenna  voltage  patterns 
of  beam  number  1 and  beam  number  2,  respectively,  measured 

with  respect  to  th?  zero  tracking  error  direction. 


FIGURE  Al:  TWO-BEAM  GEOMETRY 


From  figure  A1  it  is  seen  that 


G] (e)  - P(6  + e)  , 


62(e)  = P(-o  + e). 


(A13) 


We  now  assume  that  e is  a small  angle  so  that  P(6  + e)  and  P(-6  + e)  can 
be  approximated  by  the  first  two  terms  in  the  Mac'laurin  expansion  to 
obtain 


G^(e)  = P(<5  + e)  z P { <5 ) + eP  ’ ( 6 ) , 
G 2(e)  = P(-6  + e)  h P ( -6 ) + P 1 (— 6 ) . 


(A14) 


From,  the  symmetry  of  an  antenna  beam  it  follows  that 


P(-5)  = P(6)  , P'(-S)  - -P'(fi)  , (A15) 


so  that  (A14)  becomes 


G-j(e)  = P ( 6 + e)  = P ( 0 ) + eP'(6)  , 
G2(e)  = P(-c  r E)  = P(')  - LP'(fi). 


(A16) 
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If  we  now  substitute  (A16)  into  (A12)  the  following  expression  is 
obtained 


(A17) 


In  this  equation  Bcp  is  the  Cramer-Rao  lower  bound  on  the  variance  of  an 
unbiased  estimate  of  the  tracking  error  e when  the  conditions  illustrated 
in  figure  A1  are  satisfied. 

To  illustrate  the  type  of  result  that  is  obtained  we  shall  assume 
that  the  normalized  train  lobs  of  the  one-way  antenna  voltage  pattern  can 
be  approximated  by  a Gaussian  beam  shape  and  that  the  normalized  tv/o-way 
pattern  can  be  represented  as 


P(*) 


(A18) 


* 


l 

where  ^ is  the  half-power  beamwidth  and  K = 2 tn  2 r 1.3863.  As  will 
be  shown,  the  assumption  of  a Gaussian  beamshape  leads  to  a very  simple 
result  and  is  perfectly  adequate  for  most  radar  system  investigations. 

If  we  now  substitute  (A18)  into  (A17)  the  following  expression  is 
obtained  for  the  lower  bound 


2 

ct 

e 


> 


exp 


l 

Tiiis  means  that 


1 

2 


(A19) 
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Via  next  determine  that  value  of  6 which  minimizes  the  right  hand  side  of 
equation  (A20).  If  we  denote  this  optimal  squint  or  offset  angle  by 
6q,  then  by  differentiating  the  right  hand  side  of  (A20)  with  respect 
to  6 it  can  be  determined  that 


Furthermore 


5 (0-4247)*i 


P(6q)  *•  exp  (-1) 

p,(5-)  = ( ■ h ) exp  {-  ?>  • 


* 

If  we  now  substitute  (A2l)  into  (A20),  the  following  expression  is 
obtained  for  the  lower  bound  when  the  optimal  squint  angle  5q  is  chosen 


o > 


e - 32MR  'fin  2 


2 

JL  (.  ...e.  . ) 

J!AO  0 > * 


c2  > 
c — 


(0.1225), 


is  the  half-power  baamwidth  of  the  one-way  pattern, 

M is  the  number  of  pulses  processed  per  beam  position, 

R is  the  per  pulse  signal-to-noise  measured  relative  to  the  nose 

of  a beam. 

Most  angle  processors  come  close  to  achieving  the  lower  bound  given 
by  ( A22  ) and  for  most  analysis  purposes  it  is  adequate  to  use  this  lower 
bound  as  the  variance  of  an  estimate  of  e,  the  tracking  error  angle.  The 
next  section  will  discuss  one  technique  that  is  utilized  for  computing  the 
tracking  error  angle.  As  will  be  shown,  this  method,  called  sequential 
lobing,  asymptotically  achieves  the  variance  given  by  ( A22  ) with  increasing 
signal-to-noise  ratio,  R. 

A.  4 Amplitude  Comparison , Sequential  Lobing  Angle  Estimation 

One  method  of  obtaining  the  direction  and  the  magnitude  of  the  error 
angle  between  the  target  direction  and  the  zero  tracking  error  direction 
(which  is  usually  the  axis  of  the  antenna)  in  one  dimension  or  coordinate 
is  by  alternately  switching  the  antenna  beam  between  two  positions  (see 
figure  Al).  The  difference  in  amplitude  between  the  voltages  obtained 
in  the  two  switched  positions  is  a measure  of  the  ang"lar  displacement  of 
the  target  from  the  zero  tracking  error  direction  or  switching  axis.  The 
sign  of  the  difference  determines  the  direction  that  the  switching  axis 
must  be  moved  in  order  to  align  the  axis  with  the  direction  of  the  target. 
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Two  additional  beam  positions  are  utilized  to  obtain  the  angular 
error  in  the  orthogonal  coordinate.  Thus  a two-dimensional  sequential 
lobing  radar  might  consist  of  a cluster  of  four  feed  hours  or  ports 
illuminating  a single  antenna,  arranged  so  that  the  right-left,  up-down 
sectors  are  covered  by  successive  antenna  beam  positions.  A cluster  of 
five  feed  horns  might  also  be  employed,  with  the  central  feed  used  for 
transmission  while  the  outer  four  horns  or  ports  are  used  for  receiving. 

High  power  RF  switches  are  not  needed  in  the  latter  arrangement  since  only 
the  receiving  beams,  and  not  the  transmitting  beam,  are  switched. 

In  this  section  we  shall  derive  the  statistical  properties  of  the 

usual  sequential  lobing  angle  estimator  that  can  be  implemented  by  the  use  of 

envelope  detectors.  The  discussion  will  concern  a one-dimensional,  two-beam 
system  of  the  type  shown  in  figure  Al.  Each  beam  is  assumed  symmetric  about 
its  borcsighc  direction  and  the  beans  are  assumed  to  be  squinted  or  offset 
by  the  same  amount  on  either  side  of  the  zero  tracking  error  direction. 

The  analysis  will  first  consider  the  case  where  only  one  received  waveform 
is  processed  per  beam  position.  The  case  where  multiple  waveforms  or  pulses 

are  processed  par  beam  position  is  discussed  in  the  latter  portion  of  the 

section.  The  analysis  will  employ  the  same  notation  and  terminology  as  used 
in  the  preceeding  sections. 

A small  tracking  error  will  be  assumed  so  that  we  can  use  the  approximations 


P(6  + e)  = P(6)  + eP'(fi)  , ) 


P(-  6 + e)  = P(fi)  - eP'(<5)  , 


(A23) 
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We  shall  first  derive  an  expression  for  a. 

From  the  preceeding  sections,  it  is  seen  that  when  no  noise  is  present 
we  have 


!Z] j = AP(6  + e) 

|Z2j  = AP(-«  + e) 


(A27) 


so  that  utilizing  the  small  angle.  approximations  given  by  (A23)  we  have 
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and  thus 


a 


P(6) 

SPuT 


(A29) 


When  a Gaussian  shaped  main  lobe  is  assumed  of  the  form  discussed  in 
section  A. 4 and  when  6 is  selected  to  obtain  the  smallest  Cram/r-Rao 
lower  bound,1  then  the  constant  a has  the  value 
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a = ^ s (0.2123)*,  , 


(A30) 


where  is  the  half-power  beamwidth. 

We  shall  next  determine  the  density  function  and  the  mean  and  variance 
of  the  estimator  e given  by  (A25)  when  noise  is  present.  It  is  first  noted 
that  the  probability  density  functions  of  |Z^|  and  | | » which  will  be 
denoted  by  f^r)  and  f2(r)  respectively,  are 
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ecu  a cion  (A21). 


30 


where  0 <_ r and  IQ(*}  is  the  modified  Bessel  function  of  first  type  and 
order  zero.  Equation  (A30)  expresses  the  classical  result  that  jZ-jj 
and  j Z2 ! have  a Rice  distribution.  If  we  now  define  the  random  variable 
s,  for  0 _<  s,  with  density  function  g(s),  as 


5 * TzTT 


then  from  page  53  of  [Miller,  1954]  it  is  seen  that  g(s)  can  be  expressed 


g(s)  = -fog-g  exp  (R^2  + R^/O  + s2)J 
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where 


I.(.)  is  the  modified  Bessel  function  of  first  type  and  order  one. 


= RP2(S  + c)  , 
2 a 


A2P2(-6  + e) 

"17  ' 


RP  (-6  + e)  , 


R 


is  the  per  pulse  signal -to-noise  ratio  for  each 
beam  when  it  is  assumed  that  the  target  is  located 
at  the  boresight  of  the  beam. 


It  is  seen  that  the  random  variable  s is  the  ratio  of  two  independent  Rice 
variates  since  the  additive  Gaussian  noise  associated  with  each  pulse  or 
waveform  that  is  processed  is  assumed  independent. 

We  next  define  the  modified  signal-to-noise  ratio  RQ  as 


R = RP2(5)  , (A34) 

0 


where  6 is  the  beam  squint  or  offset  angle  zero  tracking  error  direction. 

R can  be  considered  as  the  effective  signal-to-noise  ratio  due  to  the 
o 

squinting  of  the  antenna  beams. 


l 

When  a Gaussian  shaped  main  lobe  is  assumed  and  when  o is  chosen  to  be  the 
optimal  squint  angle  defined  in  the.  previous  section,  then 

Rq  - e~ 1 P.  (n.';:.7ri}P.. 


*2 


We  now  sea  that  the  estimator  e can  be  expressed  as 


s2  - 1 

e = a “«  ) 0 < s 

s*  + 1 


where  the  density  function  of  the  random  variable  s is  given  by  (A33) 
In  order  to  find  the  probability  density  function  of  e the  following 
identities  will  be  needed. 
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where  the  latter  four  identities  utilize  small  angle  approximations  which 
are  valid  since  the  error  angle  e is  assumed  small. 

If  v/e  now  make  the  change  of  variable  from  s to  e given  by  (A35) 
and  utilize  the  above  identities  it  can  be  shown  that  the  probability 
density  function  of  e,  denoted  by  h(e),  can  be  written  as 


for 

-a  £ e 1 a • (A36) 
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Making  the  change  of  variable  c = a cos  $ it  is  seen  that  the  moment 
generating  function  can  be  expressed  as 


ir 

G(s)  = ^ (1  + Rq)  exp  (-Rq)  J sin  $ exp  (a  cos  «)I0(B  sin  $)d$ 
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(A38) 


’where 

“ = (f)R0  + as  • 

6 - \f'  -(r)2  ■ 

so  that 
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(A39) 
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The  three  integrals  given  in  equation  (A38)  are  evaluated  in  the  annex 
to  this  appendix,  and  the  moment  generating  function  G(s)  can  ultimately  be 


evaluated  as 


( 1 + Rft  ) -Rr 

G(s>  = e 0 («Y  - e-Y) 


, +■«»*;  -R0 
+ J 75 L e 


jV  + e'Y)  . 1 (eY  - e'Y,  ] , 


where  y is  given  by  (A39).  It  is  easily  checked  that  G(0)  = 1 which  must 
hold  since  h(e)  is  a density  function. 


E(e)  = 


\ds  / s=0  » 


which  after  some  algebraic  manipulation  is  evaluated  as 


1 ••  e' 


(A42) 
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Notice  tli3t  e is  a biased  estimator  for  moderate  to  small  values  of  the 
modified  signal-to-noise  ratio  Rq,  for  example, 

when  RQ  = 5,  E(£)  = (0.82 )e  , 

and 

when  Rq  = 20,  E(e)  = (0.95)e 

As  Rq  ->  » , i.e.,  as  the  noise  goes  to  zero,  it  is  seen  that  the 

estimator  e becomes  unbiased. 

One  common  way  that  is  utilized  to  remove  this  bias  is  to  use  the 
estimator  e-j  instead  of  the  estimator  e where 


]z2i2  - lzil2 

jr—  2 

\1Z\  + |Z,|  + b 


(A43 ) 


and  b is  a properly  chosen  compensation  term.  The  analysis  of  this  estimator 
and  the  selection  of  an  "optimal  value"  of  b is  quite  difficult  and  will 
not  be  pursued  further  in  this  report  [see  Lank,  Pollon,  1969]. 

Instead  we  shall  consider  the  unbiased  estimator  where 
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(A44) 
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The  estimator  has  the  disadvantage  that  an  estimate  of  RQ,  the 
modified  signal-to-noise  ratio,  is  needed  for  its  implementation.  It  is 
seen,  however,  that  the  estimator  is  not  very  sensitive  to  errors  in 
estimating  RQ  over  a large  operating  signal-to-noise  ratio  range.  From  (A40) 
it  is  seen  that  the  variance  of  e2  depends  on  the  value  of  e being  estimated. 
For  convenience  we  shall  only  compute  the  variance  of  e for  the  value 
e = 0.  However,  since  the  estimator  i is  only  suitable  for  estimating 
small  error  angles,  the  variance  for  non-zero  values  of  t will  be  closely 
approximated  by  the  variance  of  the  "s  = 0"  case. 

It  can  now  be  shown  that 


(3  + e“2s0) 
2so 
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(A45) 


Therefore,  since 


a = 


P<s) 

WTZT 


R0  = RP2(s)  , 


it  follows  that 
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Which  is  not  present 


in  the  estimator  z. 


(A46) 
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Now  from  equation  (A17)  of  the  previous  section  it  is  seen  that  is  the 

j 

Cramer-Rao  lower  bound  on  the  variance  of  an  unbiased  estimator  of  e when 
one  pulse  or  received  waveform  is  processed  per  beam  position  (i.e.,  when 
M = 1). 

Finally,  we  obtain  from  (A45)  and  (A46) 
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(A47) 


From  (A47)  it  is  seen  that  Var  (e^lc  = 0)  approaches  the  Cramer-Rao  lower 
bound  with  increasing  signal-to-noise  ratio  RQ.  However,  even  for  small 
values  of  Rq  the  estimator  is  remarkably  efficient.  For  example, 

when  R0  - 5,  Var  (e2|e  = 0)  = (1,13)B^R  . 

It  is  of  interest  to  note  that  if  the  bias  present  in  the  biased  estimator 
e given  by  expression  (A42)  can  be  "lived  with"  then  the  variance  of  this 
estimator  is  less  than  the  Cramer-Rao  lower  bound  on  the  variance  of  an 
unbiased  estimator. 


AAj  M Waveforms  or  Pulses  Processed  Per  Beam  Pos i t i on 


if"j  preceding  analy'f„  j.as  assumed  that  one  pulse  or  received  waveform 
is  pr\;:v.-s„*,d  par  beam  position.  If  M pulses  a re  processed  per  beam  position, 


then  expression  (A12)  indicates  that  the  estimator  should  be 
modified  to  read 
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(A48) 


where  |Z^|  and  | , for  j = 1 , — ,M,  are  the  envelope  detected  outputs 
from  beam  1 and  beam  2,  respectively.  This  expression,  which  gives  an  un- 
biased estimate  of  e,  has  a variance  reduced  by  a factor  of  M with  respect 
to  the  estimator  when  only  one  pulse  or  pattern  is  processed  per  beam 
p.cition. 

Frequently  the  estimator  £3  where 


(A49) 


and  k is  a normalization  constant  and  b is  a bias  compensation  term  is 
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used  instead  of  9ivsn  by  (A48).  The  reason  for  using  the  estimator 
is  to  avoid  the  possibility  of  dividing  by  a very  small  number  if  the 
amplitude  of  the  target  return  signal  is  small.  Ho  analysis  of  the 
estimator  e3  will  be  given  in  this  report.  A preliminary  analysis  is  given 
in  [Lank,  Pollon,  1969]  under  the  assumption  that  M is  large  enough  so  that 


can  be  approximated  by  Gaussian  random  variables  (utilizing  the  Central 
Limit  Theorem)  for  the  purpose  of  computing  moments.  Also,  no  analysis 
will  be  given  in  the  report  on  how  well  the  above  estimators  operate  in 
the  presence  of  a fluctuating  target  signal.  The  analysis  given  above  has 
assumed  that  the  amplitude  of  the  received  target  signal  remains  essentially 
constant  over  all  the  pulses  that  are  processed. 


E 


I 
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A. 5 Discussion 


t 


1 


$ 


Section  A. 4 of  this  report  derived  the  Cramer-Rao  lower  bour.d  on  the 
variance  of  estimating  the  error  angle  between  the  target  direction  and  the 
zero  tracking  error  direction.  A very  simple  expression  was  obtained  for 
the  case  of  a two-beam  system  where  each  beam  is  squinted  at  the  same  angle 
on  either  side  of  the  zero  tracking  error  direction.  The  expression 
obtained  assumed  that  the  main  lobe  of  each  of  the  two  beams  has  the  same 
shape  and  that  the  normalized  two-way  power  pattern  can  be  approximated  by 


* 


* 


Pfo)  = exp 

where  ^ is  the  half-power  beamwidth  and  K = 2 zn2  =.  1.3863.  The  expression 
obtained  also  assumed  that  the  two  beams  were  squinted  at  an  optimal  angle 
from  the  zero  tracking  error  direction  so  as  to  achieve  the  best  possible 
lower  bound.  It  was  shown  in  section  A. 4 that  this  optimal  squint  or  off- 
set angle  could  be  expressed  as 
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Under  the  above  conditions,  the  following  inequality  was  obtained 
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2 

(0.1226)^ 

MR 


(A52) 


where 

is  the  half-power  beamwidth  (radians)  of  the  one-way  voltage  pattern, 
M is  the  number  of  pulses  or  received  waveforms  processed  for 
each  of  the  two  beam  positions, 

R is  the  per-pulse  signal-to-noise  ratio  measured  relative 
to  the  nose  of  the  beam, 

2 

is  the  variance  of  any  unbiased  estimator  e of  the  tracking 
error  angle  e. 


The  derivation  of  equation  (A52)  assumed  that  the  noise  variates  added 
to  the  output  observables  are  independent,  zero  mean,  complex  Gaussian 
random  variables  whose  moments  are  given  by  (A3). 

Most  error  angle  processors  come  close  to  achieving  the  lower  bound 
given  by  equation  (A52)  and  for  most  analysis  purposes  it  is  adequate  to 
use  this  lower  bound  as  the  variance  of  an  estimate  of  t,  the  tracking 
error  angle.  Section  A. 5 discussed  amplitude  comparison,  sequential 
lobing  which  is  one  common  technique  used  to  accomplish  error 


angle  estimation.  Un-Jar  conditions  of  moderate  signal-to-noise  ratio, 
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i ?.  this  estimation  procedure  approaches  the  variance  given  by  the  Cramer-Rao 

j bound.  The  procedure  can  be  implemented  by  the  use  of  envelope  detectors  and 

I is  sequential  in  nature,  that  is,  the  signal  pulses  occur  sequentially  in 

time.  The  two  antenna  beams  are  produced  by  using  two  separate  horns  or 
ports  to  feed  the  antenna. 

Another  technique  that  is  commonly  used  to  accomplish  error  angle 
estimation  is  called  amplitude  comparison  monopulse.  This  procedure  is 
not  discussed  here  and  the  reader  is  referred  to  the  papers  referenced  in 
section  A.8  and  to  appendix  B.  The  method  employs  a simultaneous-lobing 
technique  in  which  the  RF  signals  received  from  two  offset  or  squinted 
antenna  beams  are  combined  so  that  the  sum  and  difference  signals  are  obtained 
simultaneously.  The  sum  and  difference  signals  are  then  processed  to  obtain 
both  the  magnitude  and  direction  of  the  error  signal.  All  the  information 
necessary  to  determine  the  angular  error  is  obtained  on  the  basis  of  a 
single  pulse,  hence  the  name  monopulse. 

An  amplitude  comparison  monopulse  system  is  less  susceptible  to 
errors  caused  by  target  cross  section  fluctuation  since  the  returns  from 
both  antenna  beams  are  received  at  the  same  time.  This  type  of  estimator 
also  achieves  a variance  which  approaches  the  Cramer-Rao  bound  with  in- 
creasing signal-to-noise  ratio.  A derivation  of  the  statistical  properties 
of  this  estimator  will  be  given  in  appendix  B. 

For  most  radar  systems  the  output  error  signal  e is  used  to  control  a 
tracking  servo  which  positions  the  two  squinted  or  offset  antenna  beams 
so  that  the  new  zero  tracking  error  direction  corresponds  to  the  predicted 
target  position  when  the  next  sequence  of  tracking  pulses  are  emitted. 
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The  direction  of  zero  tracking  error  is  changed  by  either  physically  moving 
the  antenna,  as  is  done  for  the  conventional  tracking  radar,  or  by 
electronically  moving  the  beams,  as  is  done  when  a phased  array  radar 
is  employed.  The  tracking  servo  is  usually  either  of  the  constant 
velocity  or  constant  acceleration  type.  The  analysis  of  such  a combination 
*s  given  in  [Swerling,  1954].  Of  course,  the  new  zero  tracking  error 
direction  will  not  exactly  correspond  to  the  true  target  direction  because 
of  noise  effects  in  estimating  the  previous  error  angle  (see  equation  (A52) 
and  because  of  servo  noise.  If  the  noise  power  is  large  enough  it  may 
happen  that  the  variance  of  the  estimate  of  the  tracking  error  c may  be 
so  large  that  there  is  a "significant"  probability  of  the  new  zeit)  tracking 
error  direction  being  widely  separated  from  the  actual  target  directior,. 

In  this  case  the  squinted  antenna  beams  will  not  be  positioned  properly  to 
detect  the  next  set  of  target  returns  and  the  tracking  process  will  be 
interrupted.  This  is  one  effect  that  intense  noise  jamming  attempts  to 
accomplish  in  addition  to  increasing  the  tracking  angle  variances. 

Equation  (A52)  can  be  used  to  evaluate  the  effect  of  CW  noise  jamming 
on  an  amplitude  comparison,  sequential  loting,  angle  tracker.  The  value 
of  R in  (A52)  must  be  interpreted  as  the  signal-to-jamming  power  ratio. 

The  jamming  power  is  computed  after  taking  account  of  the  beam  gain  ir. 
the  direction  of  the  jamming  emitter.  Usually  the  jamming  enters  through 
the  antenna  sidelobes  where  the  antenna  gain  is  greatly  reduced  as  compared 
to  the  main  lobe  gain.  The  signal -to- jamming  ratio  must  take  account  of 
any  ECCM  capabilities  present  in  the  radar,  for  example,  adoptive  sidelobe 
cancellation  which  places  nulls  in  the  receive  antenna  pattern  in  the 


receive  antenna  pattern  in  the  direction  of  the  noise  sources.  An 
additional  study  would  be  required  to  analyze  the  effect  of  CW  noise  jamming 
on  an  amplitude  comparison  monopulse  system  since  the  derivation  of  (A52) 
assumed  that  the  observables  contained  independent  white*  additive  Gaussian 
noise.  Such  would  not  be  the  case  in  a monopulse  system  where  the  returns 
from  both  offset  beam  positions  are  obtained  simultaneously.  Another  type 
of  jamming  that  may  be  present  is  large  amplitude,  short  duration,  pulse 
jamming.  The  effect  of  such  jamming  is  not  considered  in  this  report. 

When  two  or  more  targets  are  present,  the  radar  tracking  system  must 
be  capable  of  distinguishing  between  them  if  either  one  is  to  be  tracked 
accurately.  Without  range  or  velocity  differences,  the  conventional 
angular  tracking  methods  cannot  separate  targets  when  they  are  separated 
by  much  less  than  one  beamwidth.  The  theory  and  design  of  such  multiple 
angle  tracking  radars  is  given  in  [Lank,  Pollon,  1968]  and  [Pollon,  1968] 
and  is  not  considered  in  this  report.  These  references  derive  the  form 
of  the  data  processor  and  analyze  the  variance  of  the  multiple  angle 
estimates  and  their  relationship  to  the  Cramer-Rao  lower  bound. 
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ANNEX 


The  following  expression  is  derived  on  pages  b8-59  of  [Miller,  1964] 


f \n/Z 

/ (sin  4.) 


exp(c  cos  S1n 


e2 

777 


(n-1 )/4 


(n-U/21 


Lc  + p2 


where  Ij,( • ) is  the  modified  Bessel  function  <"f  first  type  and  order  k. 
If  we  let  n = 2 in  this  expression  and  use  the  identity 


Fit  T (1\  _ sinh(Z)  _ J_  ( ez  - e"Zl 

\2Z  h/2[l) Z 2Z  ' > • 


then  we  obtain  the  expression 


exp  (a  cos  $)  I (g  sin  $)  d;> 


ex ^ pj)-  exo  ( - V g2  + s2) 


Similarly,  if  we  let  n = 4 in  expression  [Al]  and  use  the  identity 


thsn  vve  obtain  tha  expression 
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c; 


Finally,  if  w a differentiate  equation  [A2]  with  respect  to  a 
we  obtain  the  expression 


/ 


sin  $ cos  $ exp(o  cos  +)In(8  sin  $)  d$ 


,7777  t cxp(^)+  exp  t^77) 
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o+r 


exp 


^Vo2  + 7 j - exp  |'Va2  + 62j 
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RADAR  ANGULAR  TRACKING  ERRORS:  AMPLITUDE 
COMPARISON,  MONOPULSE,  ANGLE  TRACKING 

B. 1 Introduction 

A variation  of  the  model  described  in  appendix  A for  predicting 
angular  tracking  errors  for  the  radar  tracking  the  cruise  missile  in 
the  presence  of  enemy  noise  jamming  is  described  in  this  appendix.  The 
variation  described  in  this  appendix  differs  from  that  of  appendix  A in 
what  is  assumed  about  the  radar’s  method  of  angular  tracking;  in  appendix 
A it  was  assumed  that  the  radar  uses  a method  called  "amplitude  comparison, 
sequential  lobing,"  whereas  in  this  appendix  it  is  assumed  that  the  radar 
uses  another  method  called  "amplitude  comparison,  monopulse."  This 
second  method  employs  a simultaneous-lobing  technique  in  which  RF  signals 
received  from  two  offset  or  squinted  antenna  beams  are  obtained 
simultaneously  so  that  they  have  the  same  phase.  An  advantage  of  this 
method  is  a smaller  bias  with  the  same  signal-to-noise  ratio. 

The  form  of  the  maximum  likelihood  angle  error  estimator  is  derived 
for  the  case  where  the  radar  uses  the  amplitude  comparison,  monopulse 
method.  It  is  shown  that  this  estimator  can  be  implemented  with  a phase 
detector.  In  addition,  the  first  two  moments  of  the  monopulse  maximum 

A 

likelihood  estimator  are  derived  and  a comparison  made  with  the  Cramer-Rao 
lower  bound.  The  probability  density  function  of  the  estimator  is  not 
obtained.  The  bias  of  the  monopu’se  estimator  is  less  than  the  bias  of 
the  sequential  lober  as  derived  in  appendix  A. 
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The  same  complex  signal  notation  used  in  appendix  A is  used  in  this 
appendix.  The  notation  £(•)  stands  for  the  expectation  of  the  quantity 
in  the  parenthesis  while  Re( - ) and  Iir( * ) stand  for  the  real  and  imaginary 
parts,  and  |*|  stands  for  the  absolute  value.  A superscript  p denotes 
the  transpose  conjugate  of  a vector  or  matrix  and  a superscript  * denotes 
the  conjugate  of  a complex  number. 

The  short  bibliography  started  in  appendix  A is  continued  in  Section 

B.5. 


B.2  Likelihood.  Equations  avd  Maximum  Likelihood  Estimator 

This  section  will  derive  the  likelihood  equations  and  the  form  of  the 
maximum  likelihood  estimator  when  a simultaneous-lobing  technique  is  used 
to  derive  error  angle  information.  In  this  method  the  RF  signals  received 
from  two  offset  or  squinted  antenna  beams  are  combined  so  that  the  sum  and 
difference  signals  are  obtained  simultaneously.  All  the  information 
necessary  to  determine  the  angular  error  is  obtained  on  the  basis  of  a 
single  pulse  received  through  the  two  antenna  beams,  hence  the  name 
amplitude  comparison,  monopulse  is  employed. 

It  will  be  assumed  that  each  antenna  beam  is  symmetric  about  its  bore- 
sight  direction  and  that  the  beams  are  squinted  or  offset  by  the  same  amount 
on  either  side  of  the  zero  tracking  error  direction  (which  is  usually  the 
axis  of  symmetry  of  the  antenna).  The  same  antenna  pattern  (referenced  to 
the  boresight  direction)  is  assumed  for  each  of  the  beams.  The  angle 
between  the  zero-tracking  error  direction  and  the  antenna  beam  boresight 
direction  will  be  denoted  by  <s.  A picture  of  the  geometry  involved  is 
given  in  figure  A1  of  appendix  A. 
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We  shall  first  discuss  the  case  when  one  received  waveform  or  pulse 
is  simultaneously  received  by  the  two  squinted  antenna  beams.  The  following 
two  complex  observables  are  utilized  in  the  derivation  of  the  maximum 
likelihood  estimator: 

Zi  » A[Gi  (e)  + G2  (e)]  e16  + Vp  (Bl) 

Z2  = A[Gj  (e)  - G2  (e)]  e10  * V2» 

where 

G-j ( - ) and  G,,(*)  are  the  normalized  two-way  voltage  patterns  of  beam 
number  1 and  beam  number  2,  respectively,  measured  with  respect 
to  the  zero  tracking  error  direction, 

A is  the  amplitude  of  the  received  pulse  normalized  to  the  antenna 
beam  boresight  direction, 

e is  the  angle  between  the  target  direction  and  the  zero  tracking 
error  direction, 

e is  the  phase  of  the  received  signal. 

The  Vm  are  the  signals  caused  by  system  noise  and  are  assumed  to  be  zero 

mean,  complex  Gaussian  random  variables  which  are  mutually  independent.  In 

addition,  the  in-phase  and  quadrature  (real  and  complex  parts)  of  each  Vm 

2 

are  assumed  independent  with  the  same  variance  a . Under  these  conditions, 
the  complex  Gaussian  random  variables  Vm  have  the  properties  [Reed,  1962] 

E<v0 ' **  • W//  - ■ (B*> 


and  all  other  first,  second,  third,  and  fourth  order  moments  are  zero,  where 
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2 

0 Is  the  variance  of  the  In-phase  and  quadrature  noise  variates, 

* is  the  complex  conjugate  operation. 

An  amplitude  comparison  monopulse  system  Is  less  susceptible  to  errors 
caused  by  target  cross  section  fluctuation  than  a sequential  lober  since  the 
returns  from  both  antenna  beams  are  received  at  the  same  time.  This 
technique  is  particularly  useful  where  pulse-to-pulse  amplitude  fluctuations 
due  to  target  variations  or  interference  signals  can  degrade  conical  or 
sequential  scanning  tracking  techniques.  The  phase  shifts  through  the  RF 
and  IF  portions  of  the  system  must  be  carefully  equalized  to  maintain  the 
equal  phase  relationship  in  both  channels  that  is  indicated  by  equation 
(Bl).  In  addition,  the  tolerances  between  the  receiving  horns  and  the 
comparator  section  of  the  feed  assembly  must  be  very  closely  controlled. 

This. appendix  will  not  discuss  the  degradation  in  tracking  performance  caused 
by  a phase  unbalance  in  the  two  channels. 

As  illustrated  in  sections  A. 2 and  A. 3,  the  likelihood  function  L 
associated  with  (Bl)  can  be  written  as 

L = InC  - r’-2  l«-K|2  (B3) 

where 


C is  a normalizing  constant. 


V 

Gj(e)  + G^'e) 

w = 

,72. 

II 

•» 

G1  U)  - G2 ( e ) ^ 

The  unknowns  in  equation  (B3)  are  e,  A and  e. 
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Now,  '.f  P(*)  is  the  normalized  two-way  antenna  pattern  measured  with 
respect  to  beam  boresight  and  & is  the  offset  or  squint  angle  associated 
with  each  beam,  then  as  shown  in  appendix  A the  following  relationships 
are  approximately  true  when  the  tracking  error  angle  c is  small, 

G]  (e)  = P(6  + c)  = P ( 6 ) + eP'(6), 

Gg ( e ) = P(-0  + e)  = P(fi)  - cP' (6)  . 

It  will  also  be  assumed  that  the  slope  of  the  antenna  pattern  is  almost 
constant  in  the  vicinity  of  the  offset  angle  6 and  that  we  can  use  the 
following  relations  for  a symmetric  pattern: 

P(S)  = P(-S),  P'(S)  = - P'(-S)  (B4) 

Using  the  above  relationship  it  is  seen  that  equation  (B3)  now  becomes: 

L ■ 1"C  - FJ-2  |M-K|2, 

where 


C is  a normalizing  constant. 


'zl 

' P(5  ) ' 

w = 

, K - 2A 

h 

.eP'(s). 

The  three  likelihood  equations  have  the  form: 

Hr  • <B5> 
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The  first  equation  of  (B5)  leads  to  the  expression 


e = 2AF(6j  Re(Z2e‘i9)  , (B6) 

where  Re( - ) stands  for  the  "real  part  of."  The  second  equation  of  (B5) 
leads  to  the  expression 


A ■ Re(2!e'1e)  , (B7) 

where  it  is  assumed  that  terms  involving  e2  and  e times  the  voltage  in  the 
difference  channel  can  be  neglected  compared  to  the  other  terms. 

Combining  (B6)  and  (B7)  we  obtain 


b Re  (Z2e~10) 
Re  (Zle~^ e) 

where  b is  defined  as: 


<B8) 


The  last  equation  In  (B5)  leads  to 


Im(Zje~16)  = 0,  (B9) 

where  Im{.)  denotes  the  "imaginary  part  of"  and  again  terms  involving  e 
times  the  voltage  in  the  difference  channel  are  neglected.  It  is  easily 
seen  that  (B9)  leads  to 


so  that  from  (B8)  the  maximum  likelihood  estimates  for  e,  which  will  be 
denoted  by  t , can  be  expressed  as 


e = 


bRe(Z2Zj*) 


TT 


(BIO) 


It  is  easily  seen  that  t is  consistant,  i.e.,  if  no  noise  is  present  then 


E(e)  = e 

If  we  note  that 


(Bll) 


Re[(Z2  + iZ 1 )Z1*]  = Re(Z2Z]*)  , 


(B12) 


then  it  follows  that  (B12)  can  be  written  as 


t 


b|Z2  + 1Z,|  [<*2  + «, ) VI  . 

IV  Re  LI  + 1zil  T^l  J 


(B13) 


Under  conditions  of  small  error  angle  e and  moderate  to  large  signal -to- 
noise  ratio  in  the  sum  channel  it  follows  that 


(B14) 


so  that  an  estimator  frequently  used  as  an  approximation  to  (B13)  Is 


no 


£ * bRe 


f(Z24  1Z,)  z,*] 

\\h -02,1  xiznj  • 


(B15) 


The  estimation  given  by  (B15)  is  very  easy  to  implement  with  a phase 
detector  u»ed  to  obtain  the  term  of  the  form  Re(ab*). 

The  estimator  given  by  (BIO)  was  first  discussed  in  this  form  by  [Mosca, 
1969].  The  estimator  computes  an  angle  estimate  once  per  pulse.  A simple 
stretching  or  “boxcar"  operation  followed  by  low  pass  filtering  supplies  d-c 
inputs  to  the  antenna  servo  amplifiers.  Additional  references  related  to 
this  estimator  are  given  by  [McGinn,  1966],  and  fr'nfstetter  and  Delong,  1969]. 

The  likelihood  function  (B3)  can  also  be  used  to  derive  the  Cramer- 
Rao  lower  bound  on  the  variance  of  an  unbiased  estimator  of  e,  the  tracking 
error  angle.  Following  the  procedure  given  in  appendix  A and  using  the 
small  error  angle  simplifications  given  by  (B4)  it  follows  that  the  lower 
bound  can  be  written  as 


2 l I i 

ae  - 4R  \ E BCR 


(B16) 


2 , 


where  a£  is  the  variance  of  any  unbiased  estimator  of  e and  is  the 
Cramer-Rao  lower  bound  and  the  signal -to>noise  ratio  R is  defined  as 


r = 

2o2 


(B17) 


Equation  (B16)  can  also  be  written  as 


Ill 


°c  - 4Rs 


where  R$  is  the  signal -to-noise  ratio  in  the  sum  channel,  i.e.. 


R . tf2M 

S A 


(B19) 


Equation  (B18)  indicates  the  bound  on  the  performance  of  an  unbiased 
angle  estimator  as  a function  of  the  sum  channel  signal -to-voise  ratio, 

R$,  and  a factor  b which  depends  on  the  slope  of  the  antenna  pattern 
and  the  beam  offset  or  squint  angle  6.  References  in  addition  to  those 
listed  above  which  discuss  the  variance  of  monopulse  estimators  are 
[Blackman,  1971],  [Lipman,  1971],  [Kerr,  1968],  [Sharenson,  1962], 
[Manasse,  1960],  [Nester,  1962]  and  [Urkowitz,  1964]. 


B. 3 Mean  and  Variance  of  the  Angle  Estimator 

We  shall  next  determine  the  mean  and  variance  of  the  estimator  given  by 
(BIO),  namely 


-2  Re(Z2Z]*)  . 


(B20) 


Since  and  have  the  joint  density  function  indicated  by  (B3),  the  sim- 
plist  way  to  determine  the  mean  of  (B20)  is  to  convert  to  polar  coordinate 


Z,  = r,e 


Z2  r2e 


(B21 ) 
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with  Jacobian  . in  these  new  coordinates,  the  estimator  becomes 
r2 

b^-cos  (g2-q}),  (B22) 


r9 

£ = - 


and  the  expression  for  the  mean  of  e becomes 


op  oo  2ir  2r 


E(e) 


■m 

0 O ~ ~ 


lr2eflri »r2*®i »e2>  d0id02dridr2 


where  £ is  given  by  (B22)  and  ftr^r^,  0-j , e^,)  is  the  joint  density  function 
in  polar  coordinates.  Rather  than  provide  all  details  of  this  integration 
we  shall  instead  indicate  the  key  steps. 

The  integration  with  respect  to  0^  and  Q2  can  performed  using  the 
identities 


2tt 


/ 


diKos(<H3)exp[acos(iH3)]  = 2^(0), 


2ir 


/ 


di|>sin(i/<-0)exp[cxcos(iH3)]  * 0, 


(B23) 


(B24) 


0 


where  1^ ( • ) indicates  the  modified  Bessel  function  of  first  type  and  order 
1.  Equation  (B23)  is  the  standard  definition  of  the  modified  Bessel  function 
and  (B24)  can  be  directly  integrated  after  making  the  change  of  variable 

Y - a COS  (iH3). 

The  integrations  with  respect  to  r^  and  r ^ are  performed  using  the 
identities 


oo 


o 


oo 


dxe 


-ax 


I] (yx) 


...  Y 


tr  iFiO*«fc» 


0 


(B25) 


(B26) 


where  is  the  confluent  hypergeometric  function.  A listing  of 

such  integrals  can  be  found  in  [Miller,  1964  appendix  1]. 

We  can  convert  the  hypergeometric  functions  to  more  recognizable  forms 
by  using  the  identities 


r.4 


= ex,  (B27) 

1F1  (1 ;2;x)  = ^£ex-l],  (B28) 

which  can  be  found  in  [GPO,  1964,  page  509]. 

The  final  result  after  performing  all  the  integrations  is 

E(e)  = e[!-exp(lRs)]  (B29) 

where  Rsis  the  signal-to-noise  ratio  in  the  suni  channel  given  by  (B19). 

This  equation  corresponds  to  equation  (20)  of  the  paper  by  Mosca,  which  was 
not  derived,  and  is  also  the  same  as  equation  (30)  in  a report  by  McAulay. 

The  McAuley  report  derives  this  equation  in  an  entirely  different  fashion. 

Equation  (B29)  shows  that  e is  a biased  estimator,  however,  the  bias  is 
negligible  when  the  signal-to-noise  ratio  R$  is  moderate  to  large.  Comparing 
equation  (B29)  above  with  equation  (25)  of  the  paper  by  Lank  and  Pollon,  which 
gives  the  analogous  expression  for  an  angle  tracking  radar  employing  sequential 
lobing,1  one  sees  that  the  bias  for  the  amplitude  comparison  monopulse  is 
much  reduced  over  that  of  a sequential  lober,  for  small  to  moderate 
signal-to-noise  ratios;  the  processing  of  phase  information  in  an  amplitude 
comparison  monopulse  rather  than  only  envelope  information  is  the  cause  of  this 
improvement. 


no  00  27T  27T 


E(£2)  * 


r1r2(cZ)f(r1  ^ e^de-jde^r^d^  (B30) 


0 0 0 0 


Also  derived  in  appendix  A. 
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and  is  easily  seen  to  be  infinite.  The  quickest  way  to  show  tnis  is 
to  note  that  the  existance  or  nonexist  nee  of  the  second  moment  is 
not  affected  by  the  fact  that  a signal  is  present.  Thus,  when  (B30) 
is  integrated  with  only  noise  present,  the  result  is  that  the 
second  moment  is  infinite.  The  basic  reason  is  that  the  denominator 
of  the  integrand  has  a singularity  at  r^  = 0.  This  result  was  noted 
in  the  paper  by  Mosca. 

The  effect  of  the  singularity  is  to  make  the  "tails"  of  the 
density  function  of  e slightly  too  high  so  that  the  second  moment  does  not 
exist.  A similar  phenomenon  happens,  for  example,  for  the  quotient  of 
two  normal  distributions  which  results  in  a Rayleigh  distribution  for 
which  no  moments  exist. 

However,  the  behavior  of  the  estimator  t can  be  very  closely 
approximated  near  the  value  e=0  (which  is  of  most  interest)  by  a 
random  variable  with  a finite  second  moment. 

It  can  be  shown  that,  under  most  conditions,  |Z^|  can  be  approximated 


|Z,|2  * 4A2P2(6)  . |E(Z,)|2  . 


(B31) 


Thus,  as  an  approximation  to  e we  can  utilize  the  computationally  simpler 
estimator  t given  by 


t = b 


Mz2z*) 


(B32) 


1 16 


to  evaluate  the  operation  of  f.  It  is  seen  that  all  moments  of  ? 
exist.  The  computation  of  the  first  moment  of  f results  in 

E(t)  = e,  (B33) 

where  we  utilize  the  fact  that  and  Z2  are  independent  and  that 
E(Re(Z2Zj*))  = Pe(E(Z2Zi*)).  The  variance  of  t is 


,ar(,)  * ^ [’ + i + i ] • 


(B34) 


so  that  for  small  error  angles,  e,  and  moderate  to  large  signal-to-noise 
ration  R&,  we  can  approximate  the  variance  by 


Var(t)  = — , 
4R_ 


(B35) 


which  is  the  Cramer-Rao  bound  given  by  (B18).  The  computation  of  (334) 
uses  the  fact  that 


E [te<z2z,*)]2  = | ReE(Z2Z2*Z1Z1*)  + 1 ReEtf^zy,*),  (B37) 
and  the  moment  theorem  of  Reed. 

The  key  assumption  for  the  approximation  of  e given  by  (B32)  is  that 
the  angle  estimator  is  operating  in  a situation  where  the  total  power  in 
the  sun  channel  can  be  closely  approximated  by  the  signal  power  only. 

Of  course  (B32)  cannot  be  used  in  an  operational  estimator  since  the  signal 
power  is  actually  unknown  to  the  angle  estimator  and  only  the  total  power 
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in  the  sun  channel  can  be  measured.  Operationally,  the  way  to  produce 
an  angle  estimator  for  which  all  moments  exist  (similiar,  for  example, 
the  angle  estimator  when  sequential  lobing  is  employed)  is  to  add  a 
small  positive  bias,  a,  to  the  denominator  of  (B20)  so  that  it  becomes 


Re  (I2Z}) 

(12,1  + a)2 


(B38) 


All  moments  for  exist  since  the  denominator  is  always  greater  than 
zero.  The  analysis  of  ea  is  quite  involved  and  has  not  been  undertaken. 
Depending  on  the  operating  range  of  R$  expected,  the  normalizing  constant 
b can  be  chosen  so  that  ea  is  unbiased.  The  use  of  such  a bias  term  might 
be  desirable  in  order  to  limit  the  estimated  value  of  e when  several  signals 
interfere  in  the  sum  channel. 


It  can  be  shown  (see  [Cramer,  1946],  pp  500-504)  that  as  R$  increases 
the  maximum  likelihood  estimator  e given  by  (B20)  becomes  asymptotically 
unbiased  and  is  an  asymptotically  normal  and  asymptotically  efficient 


estimator  of  e,  i.e.,  the  variance  approaches  the  Cramer-Rao  lower  bound. 

A A 

T is  means  that  e and  also  t can  be  approximated  by  a normal  distribution 

b2 

with  mean  e and  variance  BCR  = ^ in  the  vicinity  of  the  mean.  Precisely 
how  far  out  on  the  tails  of  the  normal  density  function  this  approximation 
holds  has  not  been  investigated,  but  the  approximation  should  be  very 
adequate  in  the  vicinity  of  e.  This  is  the  region  of  interest  in  most 
applications  which  involve  the  correction  of  smal1  tracking  errors. 
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APPENDIX  C 

THE  DIFFERENTIAL  EQUATIONS:  AN  ANALYTIC  APPROACH 


The  overall  mathematical  structure  of  the  single  ship/mult'iple  cruise 
missile  engagement  model  Involves  the  solution  of  a sequence  of  systems  of 
differential  equations 

m j ( t+6 1 d ) = mj (t jd)  • |l-«*z  P1  j (tjd)  kij-(t|d) J (Cl) 


for  t£  < t < t^,  k+1  < j <_  M and  all  d.  There  is  one  such  system  for 
each  of  M+l  time  intervals  [t£,  t^j).  The  assumed  ship's  weapon  assignment 
logic  led  to  expressions  for  the  firing  probabilities  of  the  form 


Pp  (*M) 


1 

m.j(t|u)  if  d=u  and  0 <.  t < t1 , 
undefined  if  d 1 u and  0 £ t < , 

0 otherwise 


and,  for  2±  j < M,  (C2) 


[l_mj-l ( t J d ) /mj ( 1 1 d ) ] 1 < i < W, 

_ (t } d)  = all  d and  0 <_  t < t. 

H j J 

0 1 <_  i < W,  all  d and  t > t . 

J 


As  has  been  explained,  equations  (Cl)  and  (C2)  can  be  used  in  tandem 
(i.e.,  by  what  has  here  been  called  the  "bootstrap  technique")  to  compute 
the  cruise  missile  survival  probabilities  m ^ ( • ( d ) . 


22 


It  Is  interesting  to  note  that  it  is  possible  to  analyze  the  systems 
(Cl)  and  (C2)  to  a certain  extent  using  only  analytical  methods.  For 
some  purposes  - - such  as  obtaining  information  about  the  structure  of  the 
solution  to  (Cl)  and  (C2)  but  not  necessarily  the  solution  itself  — these 
analytic  methods  may  suffice  and  be  more  convenient  than  the  usual  solution 
procedure  involving  numerical  bootstrapping.  The  purpose  of  this  appendix 
is  to  briefly  describe  the  analytic  approach  to  (Cl)  and  (C2). 

To  describe  the  analytical  approach  to  working  with  (Cl)  and  (C2),  it 
suffices  to  discuss  the  case  of  one  of  the  time  intervals  — say 
tj|  < t<  tj^+,  for  some  fixed  k with  1 < k<N  — and  one  of  the  damage 
states  d;  the  analysis  for  the  other  time  intervals  and  damage  states  is  the 
same. 

Writing  the  differential  equations  (Cl)  in  the  more  customary  form 

dm.(t|d)  W 

— * = -m.(t|d)*z  P.* i(t(d)-k.  -(t |d) 

dt  J i=l  1J 


and  inserting  the  expressions  (C2)  for  the  firing  probabilities  p..(t|d) 

^ J 

gives 


dm.(t|d)  W 

= [mi_1(t|d)-m.(t|d)]*z  k..(t|d) 

dt  J 1 J i=l 

= aj(t|d)-Cnij_1  (t|d)-mj(t|d)] 

for  j=k+1,  •••,  M and  t£  < t<  t^,  where 

W 

a.(t|d)  = z k..(t|d) 


(C3) 
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for  j=k+l,  •••,  M and  < t<  tfc+1.  The  initial  condition  for  the  system 
isthe  one  which  applies  to  the  time  interval  [t£,  tk+j]  namely1 

m j ( 1 1 d ) = m^t^  jd)  when  t«t£,  k+1  < j<H. 


The  equations  in  the  system  (C3)  have  a special  fomu  they  are  often  called 
"first-order  linear"  equations. 

It  is  convenient  to  write  (C3)  and  (C4)  in  a more  compact  form  by 
using  vector  notation.  Put 


*nT(t|d) 


Aft) 


(ak+1 (t|d)  -ak+1(t|d)  0 ••••  0 

0 a(c+2^tld^  “ak+2^d)  **•*  0 

00  0 aM(t|d) 


lThe  values  mj(tk|d)  may  be  computed  from  (18)  in 


the  usual  way. 
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Then  (C3)  and  (C4)  may  be  written 


~ltft)  * A(t)"m(t),  (C5) 

~m(t|d)  *-  m(tk|d)  when  t=tj^. 

If  A(t)  were  a constant  matrix  — say  A(t)=A  for  all  t in  the  range 
t£  <.  t < tk+1  — then  the  unique  solution  to  the  system  (C5)  would  be1 

m?t|d)  = exp[A{t-tk)]-"m(tk|d)  (t£  < t < tk+1). 

Since  A is  a band-diagonal  matrix,  the  computation  of  the  exponential  matrix 
would  be  straightforward. 

More  generally,  when  A(t)  is  not  a constant  matrix  (this  is  the  case  most 
likely  to  arise  in  practice),  it  may  be  shown  (set,  e.g.,  [Bellman,  1970]) 
that  the  unique  solution  to  (C5)  is 

”m(t|d)  = X(t)-m(tk|d)  (tk<t<tk+1) 

where  X(t)  is  the  unique  matrix  satisfying 

^ ■ A(t)  X(t), 

*(<)  ■ I- 

xIn  particular,  m(tk+1|d)  would  be  given  by  m(tk+1|d)  = exp[A(tk+1-tk)]^m(tk|d) 
and  these  values  of  the  rn^j^  |d)  would  be  used  in  computing  the  mj ( tk+i  I d) 
in  the  usual  way  for  the  next  time  interval. 
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The  matrix  X(t)  may  be  computed  using  the  so-called  “method  of  successive 
approximations": 

X(t)  = lim  Xn(t) 
n+-*> 


where 


x0(t)  = I 

V‘>  ' 1 A(s)Xn(s)ds 


n=0,l ,2* • • . 


Notice  finally  that  the  above  described  "analytical  approach"  does  not 
completely  avoid  the  use  of  numerical  techniques  because  the  method  of 
excessive  approximations  is  inherently  numerical  in  character  and  because 
the  computation  of  the  kernel  matrix  A(t)  (for  use  in  C6)  will  normally  be 
accomplished  numerically. 


